Calhoun 

iniutuiiaiul  Archive  Nil  vdl  Pat(<)radua(«  School 


Calhoun:  The  NPS  Institutional  Archive 
DSpace  Repository 


Theses  and  Dissertations 


Thesis  and  Dissertation  Collection 


1976 

Social  cost  of  oil  pollution. 

Mudjiardjo,  Hanny  Susmono 

Monterey,  California.  Naval  Postgraduate  School 


http://hdl.handle.net/10945/17945 


Downloaded  from  NPS  Archive:  Calhoun 


DUDLEY 

KNOX 

LIBRARY 


Calhoun  is  s  project  of  the  DudlEy  Knox  library  at  NPS,  furthering  the  precepts  and 
goals  of  open  government  and  government  transparency.  All  information  contained 
herein  has  been  approved  for  release  try  the  NPS  Public  Affairs  Officer. 

Dudley  Knox  Library  /  Naval  Postgraduate  School 
411  Dyer  Road  /  1  University  Circle 
MontereVr  California  USA  93943 


http://www.nps.8du/ljbrary 


SOCIAL  COST  OF  OIL  POLLUTION 


Hanny  Susmono  Mudjiardjo 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


SOCIAL  COST  OF  OIL  POLLUTION 
by 

HANNY  SUSMONO  MUDJIARDJO 

MARCH  1976 

Thesis  Advisor;  K.  Terasawa 


Approved  for  public  release;  distribution  unlimited. 

T174150 


Unclassitted 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whmn  D*tm  Snffd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  number 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT’S  CATALOG  NUMBER 

4.  title  5ubfa/«; 

Social  Cost  of  Oil  Pollution 

5.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Master's  Thesis 

March  1976 

A.  performing  ORG.  REPORT  NUMBER 

7.  AUTHORr»; 

HANNY  SUSMONO  MUDJIARDJO 

I.  CONTRACT  OR  GRANT  NUMBER^*; 

9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  93940 

10.  program  ELEMENT.  PROJECT,  TASK 

AREA  A  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Naval  Postgraduate  School 

Monterey,  California  93940 

12.  REPOAT  DATE 

March  1976 

IS.ggUMBER  OF  PAGES 

14.  monitoring  AGENCY  NAME  A  ADORCSSO/  dUimnmt  from  Controiling 

IS.  SECURITY  CLASS.  Ca<  If>l0  riport) 

Unci assif ied 

1$«.  DECLASSI  FI  cation/ down  grading 

SCHEDULE 

t«.  OISTRIBUTION  STaTSMEMT  (oI  £/>/•  F»port) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  ^’o/  thm  in  BioeA  20,  if  diUmfmtt  from  Rmport) 

18.  supplementary  NOTES 

19.  KEY  WORDS  (Continuo  on  fovorao  H  nocommorj  ond  Idontifjr  by  AiocA  nuaiA«0 

/■ 

20.  ABSTRACT  (Contirmo  on  r«r«r««  «id«  U  nmcoommy  md  idontity  by  block  numborj 

The  purpose  of  this  thesis  is  to  investigate  a  mehtod  of  improving  decision 
making  relative  to  the  problems  created  by  oil  spillage.  Many  countries 
around  the  world,  including  Indonesia,  are  plagued  by  increasing  pollution 
from  these  spil 1 s. 

This  thesis  uses  a  simulation  to  consider  the  spread  and  damage  caused  by 
oil  spills  using  data  from  San  Francisco  Bay.  A  projection  of  social  costs 

DD  I  jam*?!  1473  cdition  OF  I  nov«*i$o*$olete  unclassified 

(Page  1)  S/N  0102-014-  6601  | 


SECURITY  CLAillFICATlOM  OF  THIS  PAGE  Dmtm 


Unci  assified _ _ 

f^lCU^^TV  classification  of  This  Of’tm  Entmro^d 


20.  from  these  spills  has  been  made. 

Formulation  of  a  methodology  for  deriving  the  social  cost  of  oil  spills  is  a 
prerequisite  in  reaching  optimal,  rational  decisions  in  managing  oil  pollution. 
Such  decisions  may  include  the  establishment  of  a  fine  structure,  determination 
of  the  required  level  of  clean-up  and  identification  of  socially  significant 
spills. 


DD  Form  1473 
,  1  Jan  73 

S/N  0102-014-6601 


Unclassified 


security  classification  of  this  PAGEr»^*r  D*tM  Entmrmd) 


SOCIAL  COST  OF  OIL  POLLUTION 


BY 


HANNY  SUSMONC  MUDJIAROJC 
MAJOR,  INDONESIAN  ARMY 

LECTPICAL  ENGINEER,  BANDUNG  INSTITUTE  OF  TECHNOLOGY, 


SUBMITTED  IN  PARTIAL  FULFILLMENT  OF  THE 
REQUIREMENTS  FOR  THE  DEGREE  OF 

MASTER  OF  SCIENCE  IN  GPERATICNS  RESEARCH 

FROM  THE 

NAVAL  POSTGRADUATE  SCHOOL 
MARCH  1976 


1965 


I 


OUDtEY  KNOX  OBRARY 

NAVAL  POSTGRADUATE  SCt 
r/lONTEREY,  CALIFORNIA  939*^ 


ABSTRACT 


Tl-E  PURPOSE  OF  THIS  THESIS  IS  TO  INVESTIGAT 
IMPROVING  DECISION  MAKING  RELATIVE  TO  THE  PR08 
BY  OIL  SPILLAGE.  MANY  COUNTRIES  AROUND  THE  WOR 
INDONESIA,  ARE  PLAGUED  BY  INCREASING  POLLUTION 
SPILLS. 

THIS  THESIS  USES  A  SIMULATION  TO  CONSIDER  T 
DAMAGE  CAUSED  BY  OIL  SPILLS  USING  DATA  FROM  SA 
BAY.  A  PROJECTION  OF  SOCIAL  COSTS  FROM  THESE  S 
MADE. 

FORMULATION  OF  A  METHODOLOGY  FOR  DERIVING  T 
OF  OIL  SPILLS  IS  A  PREREQUISITE  IN  REACHING  OP 
NAL  DECISIONS  IN  MANAGING  OIL  POLLUTION.  SUCH 
INCLUDE  THE  ESTABLISHMENT  OF  A  FINE  STRUCTURE, 
OF  THE  REQUIRED  LEVEL  OF  CLEAN-UP  AND  IDENTIFI 
SOCIALLY  SIGNIFICANT  SPILLS. 


E  A  METHOD  OF 
LEMS  CREATED 
LD,  INCLUDING 
FROM  THESE 

HE  SPREAD  AND 
N  FRANCISCO 
PILLS  HAS  BEEN 

HE  SOCIAL  COST 
TIMAL,  RATIO- 
DECISIONS  MAY 
DETERMINATICN 
CATION  OF 
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I.INTROOUCTIGN 


f 


M4N  HAS  BEEN  POLLUTING  THE  WATERS  OF  THE  WORLD  FOR  YEARS. 
UNTIL  THE  ECOLOGY  MOVEMENT  AND  THE  PRESENT  RESOURCE  CONSERV¬ 
ATION  PROGRAMt  LITTLE  CONCERTED  EFFORT  HAS  BEEN  MADE  TO 
REDUCE  POLLUTION. 

SINCE  1970,  THE  U.S.  ENVIRONMENTAL  PROTECTION  AGENCY  HAS 
PLAYED  A  MAJOR  ROLE  IN  ATTEMPTS  TO  REDUCE  THE  FREQUENCY 
OF  OIL  AND  HAZARDOUS  SUBSTANCE  SPILLS  AND  TO  MINIMIZE  ENVI¬ 
RONMENTAL  DAMAGE  CAUSED  BY  THOSE  SPILLS  THAT  DO  OCCUR. 

IN  ADDITION,  THE  U.S.  COAST  GUARD  HAS  BEEN  INCREASING  ITS 
EFFORTS  IN  DETECTION  AND  CLEANING  OF  HAZARDOUS  SUBSTANCES 
SPILLED  INTO  BODIES  OF  WATER. 

OVER  13, COO  SPILLS  OF  OIL  (REF-B)  AND  HAZARDOUS  SUBSTANCES 
OCCUR  ANNUALLY.  SPILLED  INTO  RIVERS,  STREAMS,  COASTAL  WATERS 
ESTUARIES  AND  LAKES,  OIL  SPREAD  IN  A  MATTER  OF  MINUTES  BY 
THE  FORCE  OF  CURRENT  INDUCED  BY  THE  W IND ,S ALT N I T Y  AND  TIDES. 
SPILLS  NCT  ONLY  REPRESENT  WASTED  RESOURCES  BUT  CREATE  SOCIAL 
COSTS  TO  THE  SOCIETY  NEARBY  DIRECTLY  AND  INDIRECTLY. 

CIL  POLLUTION  IS  THE  ALMOST  INEVITABLE  CONSEQUENCE  OF  THE 
DEPENDENCE  OF  GROWING  POPULATION  ON  AN  INCREASINGLY  OIL- 
BASED  TECHNOLOGY. 

BECAUSE  OF  THE  LARGE  QUANTITIES  OFTEN  INVOLVED  IN  SPILLS, 

THE  EFFECTS  ARE  NOT  ALWAYS  COMPARABLE  TC  THOSE  CAUSED  BY  THE 
CHRONIC  POLLUTION  OF  INDUSTRIAL  AND  MUNICIPAL  DISCHARGES. 
SOME  OF  THE  EFFECTS  ARE  OBVIOUS,  SUCH  AS  POLLUTED  BEACHES, 
RIVERS  DCTTED  WITH  OIL  SLICKS,  DEAD  BIRDS  AND  FISH  . 

BUT  THE  ECOLOGICAL  EFFECTS  FROM  SPILLS  ARE  NOT  CONFINED  TO 
THE  IMMEDIATE  OR  OBVIOUS  SINCE  OVER  A  LONG  PERIOD  OIL  SPILLS 
COULD  CHANGE  THE  COMPOSITION  OF  AQUATIC  COMMUNITIES  OR  DAMA¬ 
GE  THE  ABILITY  OF  THE  SPECIES  TO  SURVIVE. 

THIS  STUDY  WILL  ADDRESS  THE  PROBLEM  OF  INCORPORATING  A  SO¬ 
CIAL  CCST  FIGURE  ON  THE  CONDITIONS  SURROUNDING  AN  OIL  SPILL. 
IT  EMPHASIZES  THE  PROBLEM  OF  GIL  POLLUTION  IN  THE  SAN  FRAN- 
CISCC  BAY  AREA  AND  ITS  RELATED  SOCIAL  COST,  A  PREDICTION  OF 
SOCIAL  CCST  HAS  SEEN  MADE  USING  COMPUTER  SIMULATION  WITH 
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PROBABLE, 


BUT  ARTIFIAL, 


INPUT  DATA. 
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II.  FREQUENCY  OF  OIL  SPILL 


A.  GIL  FCLLUTIGN  IN  THE  WATERS  OF  THE  UNITED  STATES. 


THIS  SECTION  WILL  STATISTICALLY  SHOW  GIL  POLLUTION  IN 
U.S.  WATERS:  TYPE  OF  DISCHARGE,  LCCATICN  OF  THE  DISCHARGE 
AND  THE  SOURCES  QP  THE  DISCHARGE.  DATA  PRESENTED  IN  TABLE 
ONE  AND  FIGURES  GME  THROUGH  THREE  IS  BASED  ON  THE  U.S.  COAST 
GUARC  pollution  P-iCIDENT  REPORTING  SYSTEM  (PIRSi. 

MACRO  DATA  OF  OIL  POLLUTION  TRENDS  IS  GIVE^  3Y  TABLS-l  WHICH 
CONTAINS  STATISTICS  FOR  CALENDAR  YEaRS  1971  TO  1974. 

OIL  POLLUTION  INCIDENTS  CCCURING  IN  1973  AND  1974  6Y  AREA 
AND  LGCATILN  aRE  SHOWN  IN  FIGURE-1  AND  FIGURE-2. 

THE  LEADING  SOURCES  EACH  YEAR  IN  THE  TOTAL  VOLUME  QISCHARGEO 
ARE  THOSE  WHICH  HANDLE  LARGS  VOLUME  OF  PETROLEUM  PRODUCTS, 
SUCH  AS  TANKERS,  REFINERIES  ETC. 

FlGUPE-2  SHOWS  THE  RELATIVE  VOLUME  DISCHARGED  INCLUDING 
ONLY  DISCHARGES  OVER  100,000  GALLONS  (OR  573,500  LITERS) 

IN  1973  and  1974.  i 
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POLLUTION  TRENDS  IN  ALL  U.  S.  WATERS 
CALENDAR  YEARS  1971 -1974 
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POLLUTION  INCIDENTS  BY  AREA 


GREAT 


LAKES 

% 


VOLUME 

fig-1 


Source  :  U.S.  Coast  Guard  P.I.R.S, 


(Ref-2) 


POLLUTION  INCIDFNTS  BY  LOCATION 


NUMBER 


fig-  2 

Source  ;  U.S.  Coast  Guard  P.I.R.S.  (Ref-2) 


SOURCES  OF  DISCHARGES  >100 ,000  GALLONS 


Calendar  Year  1973 


Coast  Guard  P.I.R.S.  (Ref~2) 


Source 


U.S  . 


B.  -SPILL  PROeABILITY 


Tr£  WAJCRITY  GF  SPILLS  ARE  QUITE  S.MALL  .  HOwcVER.  IT  IS 
THE  VOLUME  OF  LARGE  SPILLS  THAI  HEAVILY  IMFLUENCES  THE  SIZE 
OF  AVERAGE  SPILLAGE. 

RELIANCE  CN  THE  ANNUAL  AVERAGE  VOLUME  SPILLED  OVER  A 
PROJECTED  TIME  PERIOD  CAN  3E  QUITE  MISLEADING. 

IN  REALITY,  "HE  ENVIPONMENT AL  IMPACT  OF  OIL  SPILLS  OSPENCS 
BOTH  CN  THE  FREQUENCY  AND  SIZE  OF  SPILLS. 

IT  ALSC  CEPcNCS  ON  THE  RATE  C-F  SPILL  DISCHARGE  RELATIVE  TO 
THE  ABILITY  TO  CLEAN  UP  THE  SPILL. 

BOTH  OEV4NNEY  AND  PAULSON  (REF-2)  CONCLUDE  THAT  THE  OCCUR¬ 
ENCE  OF  A  POLLUTION  INCIDENT  IS  ESSENTIALLY  A  RANDO'^  PROCESS 
AND  CAN  BE  DESCRIBED  BY  A  POISSON  DISTRIBUTION  : 


WHERE  : 

N  =  #  OF  SPILLS 

X  =  VOLUME  HANDLED 

\  =  MEAN  SPILL  INCIDENCE  RATE  IN  H  S P  IL LS / VOL UM E 

HANDLED. 

P(N|X)  =  PROBABILITY  CF  'N*  SPILLS  OCCUR  IN3  GIVEN  X. 


THEY  CONCLUDE  THAT  : 

1.  AVERAGE  SPILL  SIZES  ARE  RATHER  MEANINGLESS  STATISTICS 
SINCE  THE  VOLUME  RANGE  IS  SO  GREAT. 

2.  THE  TRUE  IMPACT  OF  SPILLS  IS  A  FUNCTION  GF  FREQUENCY 
SIZE  AND  LOCATION. 


BRUCE  BEYAERT  (REF-5)  ESTIMATES  THE  RISK  GF  AN  OIL  SPILL  BY 
MEANS  OF  A  STATISTICAL  PROBABILITY  ANALYSIS  USING  THE  PRO¬ 
CEDURE  ILLUSTRATED  Tm  FTGURE-4, 

THIS  ILLUSTRATION  ASSUMES  THAT  AN  ADEQUATE  AND  VALID  BODY 
OF  DATA  IS  AVATLiSLE  TO  INDICATE  THE  ACTUAL  NUMBER  AND  SIZE 
OF  ACCICENTAL  GIL  SPILLS  . 

THE  CATA  WaS  USED  TQ  DETERMINE  BOTH  THE  PROBABILITY  OF  A 

THE  PROBABLE  DISTRIBUriGN  GF  SPILL  SIZE 


SPILL  EVENT  AND  ALSO 
FOR  EACH  EVENT. 

USING  THIS  ANALYSIS,  THE 
SIZE  CAN  BE  COMPUTED  FOR 
THE  RECURRENCE  INTERVAL  I 


RECURRENCE  INTERVAL  FOR  ANY  SPILL 
PANGS  OF  INTEREST. 

R  the  average  period 


BETWEEN  TWO  SPILLS  GREATER  THAN  OR  EQUAL  " 

THIS  APPROACH  IS  EASILY  UNDERSTOOD.  BUT  ITS  _  _  _ 

THE  VALID  AND  RELEVANT  OIL  SPILL  DATA,  WHICH  UNFORTUNATELY 
IS  VERY  HARD.  TO  DETAIN  IN  MOST  INSTANCES, 


OF  TIME 
4  SPECIFIED 
USE  DEPENDS 


SIZE 

CN 
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ILLUSTRATION  OF  AN  OIL  SPILL  RISK  ANALYSIS 


1.  SPILL  frequency 

USING  VALID  AND  RELEVANT  HISTORICAL  DATA,  AN  ESTIMATE  FOR 
-  N  -  IS  MACE.  THIS  ESTIMA'^ICN  IS  THE  AVERAGE  NUMBER  OF 
SPILL  EVENTS  PER  YEAR:  BASED  ON  CORRELATION  WITH  AN  APPRC- 
.  PRIATE  OPERATING  PARAMETER  SUCH  AS  THE  NUMBER  OF  TANKER  PORT 
VISITS,  THE  NUMBER  OF  OIL  TRANSFER  OPERATIONS  OR  THE  VOLUME 
OF  OIL  HANDLEC. 

2.  SPILL  SIZE  DISTRIBUTION 

USING  VALID  AND  RELEVANT  HISTORICAL  C ATA , OE T ERM INE  THE 
PROBABILITY  DISTRIBUTION  OF  ACTUAL  SPILL  SIZES. 


P*  =  PERCENT  PROBABIL I TY, THAT  A  SPILL  KILL  HAVE  A  SIZE  LESS 
THAN  OR  EQUAL  TO  S S  GIVEN  THAT  A  SPILL  HAS  CCCJRED. 

3.  RECURRENCE  INTERVAL 

FROM  THE  ABOVE, ■ R THE  RECURRENCE  INTERVAL  IN  YEARS  CAN 
dE  ESTIMATED  FOR  SPILLS  OF  VARIOUS  SIZES  BY  USING  THE  FCLLCW 
ING  RELATIONSHIP; 


R 


I 


100 

N  (  100  -  s'  ) 


FIGURE-4 


III.  THE  SPREAD  OF  OIL  SPILLAGE  IN  WATER 


A.  THE  FATE  OF  OIL 


IF  OIL  IS  SPILT  ON  LAND,  PART  OF  IT  WILL  3E  A3SOR3cO  INTO 
THE  SOIL  AND  PART  OF  IT  WILL  FLOW  OVER  THE  SURFACE  SEEKING 
A  LOW  SPOT. 

THE  FATE  OF  GIL  SPILLED  ON  THE  WATER  IS  A  VERY  COMPLEX 
SUBJECT.  MUCH  EFFORT  HAS  BEEN  MADE  TO  ESTIMATE  THE  PHYSICAL 
SPREiDiNG  AND  MOVEMENT  OF  OIL  ON  THE  SURFACE  OF  WATER  UNDER 
THE  INFLUENCE  OF  WIND,  WAVES  AND  CURRENTS. 

THE  PRECICTIVE  MODELS  DEVELOPED  SO  FAR  ARE  MOT  CAPABLE  OF 
HANDLING  SUCH  CCMPLEXlTIcS  AS  : 


1.  EVAPORATION, 01 S SOLUTION, SEDIMENT AT  ION, EMULSIFICATION, 
ALTO/PHOTO-OXIDATION  AND  3 1 CDEGR ACAT I  ON . 

2.  THF  change  IN  PHYSICAL  AND  CHEMICAL  PROPERTIES  OF 
FLOATING  OIL  RESULTING  FROM  PART  ONE. 

3.  THE  EFFECT  OF  SEA  CONDITIONS.  - 


IT  WAS  OBVIOUS  from  EXPERIMENTS  E  EXPERIENCES  THAT  SMALL 
CUANTITIES  OF  OIL,  I.E.  A  FEW  TONS  (  ONE  TON  OF  OIL  APPROX. 
EQUAL  TC  6.5  BARRELS  OR  250  GALLONS  OR  947  LITERS) 

DISAPPEAR  rapidly  FROM  THE  MARINE  ENVIRCNMENT. 

THE  GENERAL  NATURE  OF  GIL  DISAPPEARANCE  INVOLVES  SUCH  PhENC- 
MENA  AS  SPREADING,  EVAPORATION,  EMULSIFICATION,  DISSOLUTION 
I  ALTC-CXIlATIGN  and  3I0DEGRA0ATICN  (REF-10). 

IT  GRADUALLY  CISAFPcARS  THROUGH  DESTRUCTIVE  AND  DISPERSIVE 
PROCESS E5, LEAVING  an  ASPHALTIC  MASS. 

OIL  SINKS  AFTLR  ITS  DENSITY  IS  INCREASED  BY  EVAPORATION, 

BY  SCLUTICN  OF  ITS  VOLATILE  FRACTIONS, BY  INCLUSION  OF 
PARTICULATE  MATERIAL  AND  BY  OXIDATION  (203ELL  1964,  PIL?EL 
1968,  REF-5).  SINKING  MAY  HAVE  3EEN  OF  PARTICULAR  IMPORTANCE 
IN  THE  SAN  FRANCISCO  SPILL  (REF-5). 

A  GREAT  MANY  PUBLICATIONS  ON  THE  MStABCLISM  OF  HYCROCARBQNS 
IN  WATER  SUGGEST  THAT  ALL  MOLECULES  PRESENT  IN  CRUDE  OIL 
CAN  BE  ATTACKED  BY  ENZYMES  CF  THE  MI C RCGRGANI S MS .  THIS 
NATURAL  PROCESS , HOWEVER , MAY  BE  TOC  SLOW  UNDER  NORMAL 
CONDITIONS  AND  THEREFORE , UNABLE  TO  CAUSE  DESTRUCTION  OF 
LARGE  CIL  SPILLS  BEFORE  DAMAGE  IS  DONE  TO  MARINE  LIFE  AND 
COASTAL  AREAS. 

IN  THE  ALTO-aXIDATION  PROCESS,  TEMPERATURE  IS  AN  IMPORTANT 
PHYSICAL  FAC7GR.  AT  SEA,  BELOW  FIVE  DEGREES  CELCiUS,  OXIDA¬ 
TION  IS  VERY  SLOW,  IT  THUS  OCCURS  WITH  GREAT  DIFFICULTY  IN 
LATITUDES  ABOVE  75  DEGREES  NORTH  OR  BELOW  75  DEGREES  SOUTH. 
IN  ECUATGRIAL  REGIONS, THE  RATE  OF  OXIDATION  MAY  ATTAIN 
A  RATE  CF  SEVERAL  HUNDRED  GRAMS  OF  GIL  PER  CUBIC  VETER  OF 
HATER  PER  YEAR  {  G.D.  FLOODGATE MI  CRG  E I  AL  DEGRADA  .  luN'  CF 
OIL • .MAR. POLL  BUL.  3  1972,  P41-43  ) 

OIL  REMAINING  AT  SEA  FOR  THREE  MONTHS  OR  '^CRE,  LOSES  VOLUME 
continually  and  may  be  reduced  TO  asphaltic  residle  .REPRE¬ 
SENTING  AS  LITTLE  AS  15?  OF  THE  ORIGINAL  VOLUME  (SMITH  1968, 
REF-5  ). 

SIVACIEF  AND  MIKCLAj  (REF-10)  AFTER  ASSESSING  THE  ROLE  THAT 
VARIOUS  FACTORS  MIGHT  PLAY  IN  DETERMINING  THE  FATE  0=  OIL 
SPI LLS , ICENTI FI  ED  EVAPCRATICN  AS  BEING  THE  MOST  SIGNIFICANT. 
EVAPORATION  IS  ENCHANCED  BY  INCREASING  WIND  SPEED,  SEA  SUR¬ 
FACE  ROUGHNESS,  AIR  TEMPERATURE  AND  DECREASING  Oil!  FILM 
THICKNESS,  lighter  WEIGHT  OIL  EVAPORATES  FASTER  THAN  HEAVIER 
TYPES. 


I 
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8.  mcdel  of  disappearance  of  oil  spills  due  to  eva^craticn. 


TV<G  6ASIC  ASSUMPTIONS  MADE  BY  SIVADI5R  AND  MIKGLAJ  IN 
THEIR  MCDEL  : 


1  iCiRT  FPCiM  SPREADING*  THE  ONLY  PROCESS  OCCURING  TO 

-0  BIO- 

degradation  are  CONSICEKcD  NEGLIGIBLE. 


2.  OIL  CONSISTS  OF  TWO  ARBITRARY  DEFINED  PARTS  i 

A.  A  VOLATILE  FRACTION  FROM  WHICH  ALL  EVAPORATIVr 
LOSSES  OCCUR 

B.  A  RESIDUUM  FRACTION  WHICH  IS  TOTALLY  UNAFrECTED  BY 
WEATHERING. 


F  =  Cl’S'T/d  +  C2-T) 


WHERE: 


c  -  TWP  ■WFTGHT  FRACTIF>N  OF  THE  WEATHERED  OIL 
■  SAMPlI  WHICH  IS  EVAOGRATEO,  IK  =ERC::NT. 

n  r?  -  CONSTAN'^t  IN  WHICH  THE  VALUE  DEPENDS 
On'tHE  TYPE  OF  OIL,  WEATHER  aMD  WA.ER 
CONDITION. 

T  =  TIMEt  IN  MINUTES. 

AS  TIME  APPROACHES  INFINITY,  THE  VALUE  OF  'F'  CONVERGES  iQ 
20  -  22?  (REF-10). 

SIN''£  THE  REMAINING  OIL  SLICK  RESIDUE  WOULD  n^R  iPL^ 

IrAVIT^  NEARLY.THE  SA.MB  AS  SEA  ijlTES.  THtRt  I|  C...N.IO 

PORT  MECHANISM. 

C.  MCDEL  TO  DETERMINE  THE  LEEWAY  OF  OIL  oLiCKS. 

t^WAY^'^IS  DFFINih''‘As‘THE  MOVEMENrOF  OIL  SLICK  OVER  THE 
SaTER  CLE  TO  THE  ACTION  OF  THE  WIND. 

SMITH  (REF-4)  CONCLUDE  THAT  : 

ICN  : 


OSL  =  Q.0179’!‘W10  +  0.0196 
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WhEPE 


OSL  =  GIL  SLICK  LEEWAY,  IN  KNOTS. 

WIO  =  WIND  SPEED  AT  10  METERS  ELEVATION, IN  KNOTS. 


3.  THE  EXPRESSION  FOR  WIND  IN  RANGE  LESS  THAN  FIVE  KNOTS 
ShCULD  LSE: 


CSL  =  C.0199=i'W10. 


4.  OIL  SLICKS  MOVE  IN  THE  DIRECTION  OF  WIND  ACROSS  THE 
WATER  surface. 

5.  OIL  SPILL  VOLUME  WAS  NOT  FOUND  TO  AFFECT  THE  MAGNITUDE 
OF  THE  SLICK  LEEWAY, BUT  VERY  THIN  OIL  FILMS  WERE 
FOUND  TC  EXHIBIT  LITTLE  OR  NO  LEEWAY. 

6.  OIL  SLICK  LEEWAY  INCREASES  AS  A  POSITIVE  FUNCTION 

OF  SEA  STATE, BUT  THE  RELATIONSHIP  WAS  NOT  QLANTITATIVE 
LY  DEFINED. 

7.  GIL  SLICK  MOVES  ACROSS  THE  WATER  SURFACE  UNCER  THE 
INFLUENCE  OF  THE  WIND  LEAVING  A  THIN  FILM  ALONG  THEIR 
PATH. 


C.  OIL  SLICK  SPREADING  ANALYSIS. 


IT  IS  A  COMMON  OBSERVATION  THAT  OIL, WHEN  SPILLED  ON  WATER 
TENDS  TC  SPREAD  OUTWARD  ON  THE  WATER  SURFACE  IN  THE  FORM 
OF  A  THIN  LAYER.  THIS  TENDENCY  TO  SPREAD  IS  the  RESUL"  O'F 
TWO  PHYSICAL  FORCES:  THE  FQPCE  jF  GRAVITY  WHICH  CAUSES  THE 
LIGHTER  CIL  TC  SEEK  A  CONSTANT  LEVEL  BY  SPRE-AOING  H-  =  IZONTAL 
.LY,JLST  AS  IT  WOULD  ON  A  PLANE  HORIZONTAL  SOLiO  SURFACE  INC 
THE  SURFACE  TENSION  FORCE  OF  PURE  WATER, WHICH  IS  USUALLY 
GREATER  THAN  THAT  OF  THE  OIL  FILM  FLOATING  ON  WATCF. 

IN  THE  OPEN  sea, THIS  SPREADING  TENDENCY  IS  AIDED  BY  WATER 
SURFACE  MOTIONS  INDUCED  BY  WAVES, WIND  AND  TIDAL  CURRENTS. 
THESE  EXTERNAL  FORCES  MAY  BE  THE  ONLY  CAUSE  OF  SPREAD  IF 
THE  SLICK  HAS  BEEN  BROKEN  INTO  SMALL  INCEPENDEN"^  PIECES. 

THE  SPREAD  RESULTING  FROM  THIS  RANDOM  MOTION  Or  THE  SEA 
SURFACE  IS  VERY  DIFFICULT  TO  ESTIMATE.  MOST  OILS  SPILLED 
ARE  MIXTURES  CF  COMPONENTS  HAVING  VARYING  VAPOR  PRESSURE 
AND  SOLUBILITY  IN  WATER.  THE  LIGHTER  MOLECULAR  WEIGHT  ARE 
MORE  VOLATILE  AND  SOLUBLE, AND  THE  HEAVI ER/RES I OUh  IS  DENSER 
AND  MORE  VISCOUS. 

FAY  (REF-6,10)  HAS  DEVELOPED  AN  ANALYSIS  FOR  THE  SPREADING 
OF  A  ONE  DIMENSIONAL  AND  AXISYMETRIC  OIL  SLICK  AS  A  FUNCTION 
OF  TIME.  THE  ANALYSIS  IS  RESTRICTED  TO  A  FIXED  AMOUNT  OF 
CIL  IN  THE  initial  SPILL  AND  TH^  WATER  IS  FREE  OF  MOTION 
INDUCEC  EY  WIND, WAVES  AND  TIDAL  CURRENTS. 

THE  SPREADING  PROCESS  PROCEEDS  THROUGH  THREE  STAGES  IN  WHICH 
VARIOUS  FORCE  PAIRS  ARE  IMPORTANT  AND  THE  SLICK  FINALLY 
REACHES  A  TERMINAL  SIZE.  FOR  THE  AXI-SYMETRIC  CASE,  THE  =A- 
DIUS  CF  THE  SLICK  'R«  IS  RELATED  TO  TIME  AFTER  THE  SPILL 
BEGINS , 'T* , BY  THE  FOLLOWING  EQUATIONS: 
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GRAVITY-INERTIA  STAGE; 


Rj  =  (A-GVT^)  C, 


GRAVITY-VISCOUS  STAGE: 


R2=  ( 


AG 

— )  .  Q 


V 


to 


SURFACE  TENSICN-ViSCGUS  STAGE: 


1/4 


WHERE : 

Cl 

C2 

C2 

G 

A 

ao 

Po 

V 

a 

'"to 

T 

R 


1.14 

1.45 

2.30 

GRAVITY,  GR/CM**2 
^i)”  ^  Ri) 

WATER  DENSITY,  GR/CM^^J'S 
OIL  DENSITY,  GR/CN*»3 
VOLUMEt  CM=i'*3 

SPREADING  COEFFICIENT,  OYNS/CM 
WATER  KI NE NAT  I C  VISCOSITY,  CM^*2/S5C. 
TIME,  SECOND. 

RADIUS,  CM 


THE  FINAL  RADIUS  CF  THE  SPILL  : 


R 


n 


1/2 


WHERE: 

RF  =  RACIUS,  METER. 
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V  =  VOLUME, 


THE  TIME  AT  WHICH  THE  TRANSITION  FROM  THE  INERTIA  STAGE  TO 
THE  VISCOUS  STAGE  OCCURS ,• T 12 ' ,CAN  BE  FOUND  BY  EQUATING  FROM 
THE  SPILL  RADII  FROM  EQUATION  ’Rl»  AND  'R2'. 


.  /  ^2  ■ 

)  {  ^ 

■  V  r  . 

'  ^  A  G 

V  ^ 

(Jj 

SIMILARLY 

FOR  THE 

VISCOUS- 

-SURFACE  TENSICN  TRANSITION  TIME: 

IN  THIS  STUDY  THE  FOLLOWING  VALUES  ARE  USED  : 
G  =  98C  CM/SEC-  . 

O'  =  10  OYNE/CM 
P  =1.0 

(jj 

pQ  =  0.95 

LATITUDE  =  38  DEGREE  NORTH 
=  0.01  CM^  ,/SEC. 


C.  MOVEMENT  OF  SPILLED  GIL  IN  SAN  FRANCISCO  BAY  AREA. 


THE  INTENT  OF  THIS  STUDY  IS  TO  PRESENT  DESCRIPTIONS  OF 
GENERAL  SURFACE  WATER  MOVEMENTS  IN  THE  SAN  FRANCISCO  BAY 
AREA  ANC  ITS  IMPACT  ON  THE  DISTRIBUTION  AND  DISPERSAL  OF  OIL 
SPILLS  IN  THE  BAY. 

A  SIMULATION  METHOD  IS  USED  TO  PREDICT  THE  DISPERSAL  OF 
SPILLED  OIL. 

THE  SAN  FRANCISCO  BAY  SYSTEM  CAN  BE  DIVIDED  INTO  A  N^RTH  BAY 
-  CENTRAL  BAY  AND  A  SOUTH  BAY. 

PERRENIAL  ESTUARINE  CIRCULATION  IN  THE  NORTH  ANC  CENTRAL  BAY 
ANC  THE  ADJACENT  OCEAN  IS  MAINTAINED  BY  THE  SACRAMENTO  -  SAN 
JOAQUIN  RIVER  RUN  OFF.  ANNUALLY  THE  AVERAGE  SURFACE  SEAWARD 
DRIFT  EXCEEDS  FIVE  KM  PER  DAY  (REF-5). 

WIND  IS  A  MAJCR  FACTOR  IN  DETERMINING  CIRCULATION  IN  THE 
SOUTH  BAY.  LOCAL  WINDS  ARE  TOPOGRAPHICALLY  CONTROLLED  BY 
PREOCMINANTLY  NORTHWEST  WINDS  DURING  SUMMER  AND  THE  SOUTH 
EAST  WINDS  DURING  WINTER  STORMS,  WITH  THE  PREVAILING  SUMMER 
WIND  GENERALLY  TWICE  AS  STRONG  AS  THAT  IN  WINTER  (REF-5). 

THE  SOUTH  BAY  CIRCULATION  IS  RELATIVELY  SlUGGISH.THE  ANNUAL 
AVERAGE  SPEED  OF  SURFACE  OR  I  FT , REGAP OL E SS  OF  DIRECTION  IS 
BETWEEN  CNE  iiND  TWO  KM  PER  DAY.  THE  MOVEMENTS  OF  SURFACE 
DRIFTERS  IN  SAN  FRANCISCO  BAY  ARE  SHOWN  IN  FIGURc-5  (REF-5). 
AN  EXAMPLE  OF  THE  SPREAD  AND  THE  MOVEMENT  OF  OIL  SLICKS 
PLOTTED  BY  THE  COMPUTER  SIMULATION  IS  SHOWN  IN  FIGURE-6. 
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RNRLV3IS  OF  OIL  SPILL  MQUEMENT 
IN  SRN  FRRNCISCO  BRV  RRER 
fiS  PREDICTED  BV  ESTURRINE  NON  TIDRL  DRIFT 


IsEST  SOUTHERN 


IV.  FACTORS  AFFECTED  BY  THE  OIL  POLLUTION 


CRUDE  OIL  AND  OIL  PRODUCTS  SPILLED  IN  NATURE  ARE  ALTERED 
BY  cVAPCRATION,  BY  DISSOLUTION,  BY  BACTERIAL  AND  CHEMICAL 
ATTACK. 

IN  SPITE  OF  COMPLEX  PROCESSES  OCCURING  CURING  WEATHERING, 
MANY  COMPOSITIONAL  PARAMETERS  ARE  RELATIVELY  STABLE  AND  ARE 
NOT  OBLITERATED  UNTIL  AN  ADVANCED  STAGE  OF  DEGRADATION  HAS 
BEEN  REACHED  (REF-1). 

THE  STABLE  PARAMETERS  MAY  AID  THE  IDENTIFICATION  OF  AN  OIL 
POLLLTANT  AND  IN  THE  CORRELATION  WITH  ITS  SOURCE  FOR  MANY 
WEEKS  AFTER  THE  SPILL  {REF-1,2). 

THE  EFFECTS  OF  OIL  POLLUTION  CAN  BE  DIVIDED  INTO  TWO  GROUPS, 
MECHANICAL  DAMAGE  AND  TOXIC  EFFECTS. 


A,  MECHANICAL  DAMAGE. 


A  MOST  OBVIOUS  EFFECT  OF  AN  OIL  SPILL  IS  THE  AROUSED 
CITIZEN  INTEREST  WHEN  NEARBY  RESOURCES  ARE  POLLUTED  BY  THE 
SPILL.  THE  VISIBLE  RESULT  MAY  BE  AN  OILED  BEACH  WHICH  HAS 
BECOME  LESS  ATTRACT  I VE , C I SCOURAGI NG  SWIMMING  AND  FISHING. 
ALTHOUGH  OIL  COATED  BEACHES  CAN  BE  CLEANED  COSMETICALLY, 

THEY  CANNOT  BE  DISINFECTED  <R£F-2). 

IF  LARGE  CUANTITIES  ARE  SP I LL ED , TH ES E  EFFECTS  CAN  BE  IMPORT¬ 
ANT.  SEA  BIRDS  ARE  PERHAPS  THE  MOST  OBVIOUS  SUFFERERS. 

THEY  DIE  WHEN  OIL  DESTROYS  THE  NATURAL  INSULATING  QUALITIES 
OF  THEIR  FEATHERS  (REF-8)  OR  MAKES  FLIGHT  IMPOSSIBLE,  RESULT 
ING  IN  DEATH  BY  STARVATION.  DUCKS  WERE  FOUND  TO  PREEN  ABOUT 
HALF  THE  POLLUTING  OIL  FROM  THEIR  FEATHERS  WITHIN  A  WEEK, 
MOSTLY  ON  THE  FIRST  DAY  OF  OILING  .  LATER  EXPERIMENTS 
(  HARTUNG  R.  AND  HUNT,  G.S.J  WILDE.  MGMT,  1966,  P564-570) 
SHOWED  THAT  THE  INGESTION  OF  ONE  THIRD  TO  HALF  THE  AVERAGE 
QUANTITY  OF  OIL  EXTRACTED  FROM  THE  PLUMAGE  OF  BIRDS  FOUND 
DEAD  FROM  MODERATE  OILING,  PRODUCES  SERIOUS  INTERNAL  EFFECT- 
S.  BUT  THE  IMMEDIATE  CAUSE  OF  DEATH  IN  A  BIRO  OILED  AT  SEA 
IS  EXPCSLRE  OR  DROWNING. 

IN  ADDITION,  OIL  DISRUPTS  THE  BIRDS'  NESTING  GROUNDS. 

ASIDE  FROM  BIRDS,  SIGNIFICANT  IMPACT  DOES  NOT  GENERALLY 
OCCUR  UNTIL  OIL  REACHES  THE  INTERTIDAL  ZONE. 

THE  MOST  ecologically  IMPORTANT  IMPACT  COULD  OCCUR  ON  THE 
SHORE  OF  MARSHES  AND  ESTUARIES  WHICH  HAVE  EXTREMELY  HIGH 
BIOLOGICAL  PRODUCTIVITY  (REF-1,2,7). 

SMOTHERING  OF  SHORE  CRABS  AND  SEiSILE  INVERTI BRAT ES  (NON 
COMMERCIAL  WATER  ANIMALS)  CAN  OCCUR  IN  THE  UPPER  INTERTIDAL 
ZONE  OF  THE  OPEN  COAST, BUT  ARE  OF  CONSIDERABLY  LESS  ECOLO¬ 
GICAL  SIGNIFICANCE  (BRUCE  SEYAERT,  ANALYSIS  OF  OIL  ACCIDENTS 
FOR  ENVIRONMENTAL  IMPACT  STATEMENT,  REF-1  P  43). 

ON  THE  SHORE  ITSELF  MOST  SEAWEEDS  HAVE  AN  OUTER  LAYER  TO 
WHICH  OIL  DOES  NOT  CLING,  SO  MODERATE  QUANTITIES  CF  OIL  CAN¬ 
NOT  DC  MUCH  CAMAGE.  INDIRECTLY  IT  CAUSES  DAMAGE  TO  OTHER 
SPECIES  WHICH  DEPEND  ON  THESE  PLANTS  AS  FOOD. 

A  HEAVY  SPILL  WILL  BLANKET  EVERYTHING  ON  THE  SHORE  AND  CLOG 
THE  GRASSES  AND  REEDS  OF  S ALT-M ARSHcS , P ART  I CUL ARL Y  IF  THE 
OIL  HAS  BECOME  EMULSIFIED.  A  HEAVY  RESIDUE  CF  CRUDE  OIL 
REDUCED  THE  POPULATION  OF  WINKLES  (LARGE  SEA  SNAILS) 

ALTHOUGH  THE  RESIDUE  WAS  VIRTUALLY  NON  TOXIC. 

A  COATING  OF  OIL, EITHER  ON  THE  SURFACE  OF  WAT=R  OR  CN  AN 
INDIVIDUAL  PLANT, ALSO  INTERFERES  WITH  LIGHT  PENETRATION  AND 
THUS  PHOTO-SYNTHESIS  (REF-7). 


B.  TOXIC  EFFECTS. 


THE  PHYTOTCXIC  EFFECT  OF  HYDROCARBONS  HAD  BEEN  STUDIED 
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IN  MOST  DETAIL  ON  TERRESTRIAL  PLANTS  .  IT  WAS  FOUND  THAT 
TOXICITY  IN  SMALLER  MOLECULES  WAS  GREATER  THAN  IN  LARGER 
ONES. 

FISH  USUALLY  KEEP  WELL  CLEAR  OF  AN  OIL  SPILL  IF  THEY  CAN. 
SHELLFISH  ARE  AFFECTED  BY  OIL  DUE  TO  THEIR  LACK  OF  MOBILITY. 
FIVE  TO  TEN  PERCENT  CRUDE  OIL  SLOWS  THE  PUMPING  RATE  OF 
OYSTERS  THUS  AFFECTING  THEIR  FEEDING,  RESPIRATION  AND  GENE¬ 
RAL  CONDITION  (REF-7). 

THIS  HAS  MEANT  FINANCIAL  LOSSES  FOR  FISHERMEN  AND  PROCESSORS 
.  ALSC,  SOME  COMMERCIAL  SPECIES  CAN  ACCUMULATE  POTENTIALLY 
CARCENOGENIC  SUBSTANCES,  DAMAGING  THE  ORGANISM  ITSELF  OR 
MAKING  IT  UNFIT  FOR  CONSUMPTION  BY  MAN  AND  OTHER  ANIMALS. 
AQUATIC  LIFE  UNDER  CONDITIONS  OF  LONG  TERM  OR  CONTINUOUS 
EXPOSURE  TO  OIL  SPILLS  OEVcLOPES  SUBTLE  CHANGES  IN  THE  BEHA¬ 
VIOR  PATTERN  SUCH  AS  LOSING  THEIR  ABILITY  TO  SECURE  FOOD, 
AVOID  INJURY,  ESCAPE  FROM  ENEMIES,  CHOOSE  A  HABITAT,  RECOG¬ 
NIZE  TERRITORY,  MIGRATE,  COMMUNICATE  AND  REPRODUCE  (REF-S)  . 


C.  IMPACT  OF  A  DISCHARGE. 


AN  OIL  SPILL  CAN  HAVE  ADVERSE  SOC lA L, ECONOM IC  AND  EN¬ 
VIRONMENTAL  IMPACT  .  THE  SEVERITY  OF  THE  IMPACT  DEPENDS 
ON  THE  ENVIRONMENTAL  SETTING  OF  THE  AFFECTED  AREA, THE  TYPE 
AND  AMOUNT  OF  OIL  SPILT,  AND  THE  MITIGATION  MEASURES  EMPLOY¬ 
ED. 

SOCIO-ECONOMIC  IMPACTS  COULD  INCLUDE  ADVERSE  PUBLICITY  AND 
PUBLIC  CCNCERN  WHICH  COULD  TEMPORARILY  DETER  COASTAL  RECREA¬ 
TION  AND  TOURISM. 

PERSONAL  PROPERTIES  SUCH  AS  BOATS  AND  FISHING  GEAR  COULD  BE 
FOULED  BY  OIL. 

COMMERCIAL  AND  SPORT  FISHING  ACTIVITY  COULD  BE  TEMPORARILY 
STOPPED  OR  CURTAILED  IN  THE  AREA. 

THE  POSSIBLE  IMPACTS  OF  A  DISCHARGE  CAN  INCLUDE  (R5F-2)  : 


A.  HAZARD  TO  HUMANS  THROUGH  EATING  CONTAMINATED  SEAFOOD, 

B.  DECREASE  OF  FISHERY  RESOURCES,  AND  DAMAGE  TC  WILDLIFE 
SUCH  AS  SEA  BIRDS  AND  MARINE  MAMMALS, 

C.  DECREASE  OF  AESTHETIC  VALUE  DUE  TO  UNSIGHTLY  SLICKS 
OR  OILED  BEACH,  DECREASING  THE  VALUE  OF  PRIVATE  PRO¬ 
PERTIES  AND  RECREATIONAL  ACTIVITIES  AND  TOURISM, 

0.  DECREASE  IN  DIVERSITY  AND  PRODUCTIVITY  OF  SPECIES  IN 
THE  POLLUTED  AREA. 

e.  MODIFICATION  OF  H A  BIT  ATS , DEL AY  I NG  OR  PREVENTING  RE- 
CGLCNIZATION. 
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V.  THE  ECCNQi^ICS  OF  OIL  SPILLS 


I\  M/iRCH  1967  THE  OIL  TANKER  TGRREY  CANYON  FOUNDERED  OFF 
THE  SOUTHERN  COAST  OF  ENGLAND,  SPILLING  119,000  TCNS  OF 
CRUDE  OIL.  THE  OIL  SLICK  QUICKLY  SPREAD  ACROSS  NEARBY  WATERS 
ANC  FOULED  LARGE  AREA  OF  ADJOINING  ENGLISH  AND  FRENCH 
COASTS.  THE  BRITISH  GOVERNMENT  ALONE  SPENT  SB  MILLION  (REF- 
12)  CN  CLEAN  UP.  THAT  WAS  ONLY  A  PORTION  OF  TOTAL  CLEAN  UP 
COSTS.  IN  ADDITION,  THERE  WAS  EXTENSIVE  LOSS  OF  MARINE  LIFE 
AND  FOULING  OF  BEACHES  AND  COASTLINES. 


SINCE  OIL  IS  NOT  COMPLETELY  BIODEGRADABLE  OR  DOES  NOT  DETER¬ 
IORATE  RAPIDLY,  SLICKS  AND  GLOBULES  OF  OIL  ARE  VISIBLE 
THRQLGHOUT  THE  HIGH  SEA  OF  THE  WORLD.  THE  EXACT  BIOLOGICAL 
CONSEQUENCES  ARE  STILL  UNDETERMINED. 


IT  IS  TE 
SHOULD  6 
COME  IF 
THE  BENE 
FROM  THI 
INDEED  8 
IT  IS  CR 
GET E PM  IN 
PRODLCTI 
OF  OIL  S 
SECTION 


MPTING  FOR 
E  STOPPED, 
THE  COSTS 
FITS. 

S  POINT  OF 
£  SOCIALLY 
ITICAL  TO 
E  THE  SOCI 
VE  INQUIRE 
PILLAGES. 
ARE  DEVOTE 


PEOPLE  TO  ASSERT  THAT  ALL  POLLUTION 
BUT  THE  SOCIETY  WILL  HAVE  LESS  REAL  IN- 
OF  TOTAL  ELIMINATION  OF  POLLUTION  EXCEED 

VIEW  SOME  LEVEL  OF  OIL  POLLUTION  MAY 
DESIRABLE. 

DEVELOP  A  PROCEDURE  AND  METHODOLOGY  TO 
AL  COST  OF  AN  OIL  SPILL  SO  AS  TC  CONDUCT 
S  IN  DERIVING  THE  SOCIALLY  OPTIMAL  LEVEL 
THIS  SECTION  AS  WELL  AS  THE  FCLLOWING 
D  TO  SUCH  DEVELOPMENT. 


IF  A  SPILL  OCCURS,  THE  DIRECT  LOSS  OF  A  PRODUCT  TO  THE  ECO¬ 
NOMY  MAY  BE  MEASURED  IN  TERMS  OF  THE  MARKET  VALUE  OF  THE 
PRC  CLCT  . 

THE^INDIPECT  COST  ASSOCIATED  WITH  THE  ENVIRONMENTAL  DAMAGE 
IS  MUCH  MORE  COMPLICATED. 

THE  INCIFECT  COST  IS  THE  FUNCTION  OF  SIZE,  FREQUENCY, 
LOCATION  AND  TYPE  OF  OIL  SPILLED  INTO  THE  WATERS. 

THE  SOCIAL  COST  IS  DEFINED  AS  VALUATION  OF  LOSSES  IN  REAL 
GOODS  ANC  SERVICES  RESULTING  FROM  THE  CIL  SPILL. 


IN  THE  ABSENCE  OF  ANY  CLEAN  UP  PROCEDURES  (DETECTION, 
CHEMICALS  EQUIPMENT  ETC),  THE  SOCIAL  CCST  OF  A  SPILL  COULD 
BE  DEFINED  AS  THE  SUM  OF  DIRECT  AND  INDIRECT  COSTS. 


C(S)  =  A(S)  +  B(S) 


WHERE: 

A,  THE  DIRECT  COST  IS  A  FUNCTION  OF  SPILL  SIZE  AND 

e,  THE  INDIRECT  COST  IS  THE  FUNCTION  OF  SIZE, 
FREQUENCY  ,  LOCATION  AND  TYPE  CF  OIL. 

TYPE  CF  OIL  CHARACTERIZED  THE  TOXICITY  TO  MARINE 
LIFE. 


A.  GENERAL  ECONOMIC  ANALYSIS 
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ASSUME  that  there  ARE  N  COMMOOTTIES  IN  THE  ECONOMY  WHICH 
CAN  SERVE  AS  GGGOS  WHERE  THE  COMMOOITItS  APE  DEFINED  FOR  A 
PARTICULAR  DATE  AND  PLACE  30  THAT  A  SINGLE  PHYSICAL  COMMO¬ 
DITY  OELIVEREC  AT  TWO  DIFFERENT  DATES  CR  TWO  DIFFERENT 
PLACES  WOULD  BE  CONSIDERED  DIFFERENT  ECONOMIC  CCM '^GOI T  I  E  S . 
ASSUME  N  IS  FINITE,  AND  THE  QUANTITIES  OF  ANY  CGMVOCirY  ARE 
assumed  PERFECTLY  DIVISIELE. 

A  PARTICULAR  BUNDLE  OF  COMMODITIES  IS  SUMMARIZED  BY  THE 
COLUMN  VECTOR  *X': 


X=  (  Xj^tX2»»»»»Xj^  )' 

THIS  VECTOR  IS  DEFINED  ON  EUCLIDEAN  N-SPACE,-  i  REFEREC 
TO  AS  CCMMCDITY  SPACE. 

PRICE  IN  THE  ECONOMY  ARE  SUMMARIZED  BY  ROW  VECTOR  P  : 


P  —  (  P  ^  2  ^ 


PRICE  are  non  (NEGATIVE  AND  AT  LEAST  ONE  PRICE  IS  NON  ZERO. 
THE  PRICE  CAN  BE  NORMALIZED  AND  ONE  POSSIBLE  N OP M AL I Z AT  I  ON 
IS  THAT  OF  MEASURING  PRICES  SO  THAT  THEY  SUM  TO  UNITY. 


N 


.2 


5- 


1 


EACH  CF  THE  FIRMS  IN  THE  ECONOMY  MUST  SELECT  LEVELS  OF 
INPUTS  AND  OUTPUTS  SUBJECT  TO  THE  AVAILABLE  TECHNCLCGY,  SO 
AS  TG  MAXIMIZE  PROFITS.  FOR  EXAMPLE  THE  FIRM  f  MAY  CHOOSE 

INPUT-OUTPUT  VECTOR  Y^  IN  THE  COMMODITY  SPACE  : 


THE  FRGCLCTICM  POSSIBILITIES  SET  A  SUBSET  CF  COMMODITY 

SPACE  IN  which: 


f  =1,^ 


IT  IS  ASSUMED  THAT  £aCH  PRODUCTION  POSS  I SIL I T I ES  SET  IS  IN¬ 
DEPENDENT  OF  THE  INPUT/OUTPUT  VECTOR  CHOSEN  BY  CTHER  FIRMS 
AND  OF  THE  CONSUMPTION  CHOICES  OF  CONSUMERS. 

THE  ECONCMY  -WIDE  I-O  VECTOR  Y  ,  IS  OBTAINED  BY  SUMMING 
ALL  INDIVIDUAL  FIRM  I-O  VECTORS: 


Y 


F  f 
2Y 
{•I  ^ 


) 


BY  SUMMATION,  INTERMEDIATE  GOODS  CANCEL  OUT,  SO  ONLY  FINAL 
QUTPLTS  (MEASURED  AS  POSITIVE)  AND  PRIMARY  RESOURCES  (MEA¬ 
SURED  AS  NEGATIVE)  APPEAR  IN  Y  , 

THE  ECONCMY-WIOE  PRODUCTION  POSSIBILITIES  SET  T  IS  OBTAINED 
BY  SLMMING  ALL  FIRM  PRODUCTION  PGSS I B  I  L IT  I ES  SETS  : 


F_  £ 

Y  €  T  =  IT 

f«t 
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ASSU^PTICNS  : 

1.  T  IS  CONVEX:  Y,Z  - >  ctY  +  (1-  a)Z  et  0  <  a<  1 

2,  IT  IS  IMPOSSIBLE  TO  PRODUCE  OUTPUTS  USING  NC  INPUTS. 

2.  OUTPUT  AND  INPUT  CANNOT  BE  REVERSED. 

4.  IT  IS  POSSIBLE  TO  USE  INPUTS  ANO  PRODUCE  NO  OUTPUT. 
INPUTS  BEING  FREELY  DISPOSABLE. 

SINCE  OUTPUTS  ARE  MEASURED  AS  POSITIVE  AND  INPUTS  AS 
NEGATIVE,  THE  PROFIT  OF  FIRM  f  : 

f  -  -f  f 

n  =  p.Y^  =  2P.  .Y; 

J-t  J  J 


TOTAL  PROFIT  n  IS  MAXIMIZED  WITHIN  f  IFF  ALL  FIRMS  MAXIMIZE 
THEIR  INDIVIDUAL  PROFITS  WITHIN  THEIR  PRODUCTION  POSSI¬ 

BILITIES  SETS  . 

EACH  OF  THE  CONSUMER  H  IN  THE  ECONOMY  MUST  SELECT  LEVELS 
OF  PURCHASES  SUBJECT  TO  A  BUDGET  CONSTRAINT. 

CONSUMER  h  SELECTS  A  CONSUMPTION  VECTOR  C^  : 


THE  TASTE  OF  CONSUMERS  ARE  SUMMARIZED  BY  THE  PREFERENCE 
RELATION,  ASSUMED  CONVEX  ANO  CONTINUOUS.  ALSO  ASSUMED  THAT 
THE  PREFERENCE  RELATION  FOR  ANY  CCNSUMERS  IS  INDEPENDENT  OF 
THE  CCNSUMPTICN  CHOICES  OF  CTHER  CONSUMERS. 

THE  EUCGET  CONSTRAINT  3  : 


-  V.  ^  u 

B  =  P.C"  =  2P  .C^ 
J=i  j  j 


total  consumption  levels  for  the  ECONOMY  C  ,IS  OBTAINED  8Y 
SUMMING  ALL  INDIVIDUAL  CONSUMER  CONSUMPTION  VECTORS  : 


Kii 


IF  A  REPRESENT  TOTAL  RESOURCES  FOR  THE  ECONOMY,  THEN 
WEALTH  OF  THE  ECONOMY  W  : 


W  =  P.A 


A  COMPETITIVE  EQUILIBRIUM  IS  DEFINED  AS  A  SITUATION  IN  WHICH 
PRICE  VECTOR  SATISFIES  : 
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L 


<t.(p*)  <  0 


i  =  1,2,  N 


(non  positive  excess  demand  function) 


P'l  *  4>i  (  P*)  =  0 


t  =  1,2,. .  .  ,N 


THE  PROFIT  MAXIMIZING  I -0  VECT0=^  OF  EACH  FIRM  IS  SUMMARIZED 

_*  _1*  '  _2*  _F* 

Y  =:  (  Y  ,Y  ,  .  .  .  ) 

THE  EQUILIBRIUM  CONSUMPTION  VECTOR  OF  EACH  CONSUMER  : 

_*  — H* 

C  =  (  C  ,C  ,  .  .  .  ,C  ) 

PR’CFIT  maximizing  subject  to  THE  AVAILABLE  TECHNOLOGY  AND 
POLLLTICN  LEVEL  : 


_*  _f*  _*  _f 

=  P  .Y  >  P  .Y 


FOR  ALL  Y^  €  T 


A  PARETO  .GPTi MUM  IS  A  SET  OF  CONSUMPTICN  VECTORS  : 


1*  2* 

(  C  » C  » 


) 


WHICH.  IS  CONSISTENT  WITH  THE  TECHNOLOGY  AND  BUDGET 
AND  FOR  WHICH  THERE  EXISTS  NO  OTHER  SET  CF  CONSUMPTION 
VECTORS,  SUCH  THAT  NO  CONSUMER  IS  WORSE  OFF  AND  AT  lEAST 
ONE  IS  SETTER  OFF. 

IN  THE  CASE  OF  ONE  CONSUMER  AND  ONE  PRODUCER  IS  ILLUSTRATED 
IN  FIG-  7. 


FIG-7 
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BY  ASSU^PTICM  THE  PREFERENCE  SETS  AND  THE  PROOUCTIGN  POSSIEI 
LITIES  ARE  CONVEX,  COMPETITIVE  EQUILIBRIUM  IS  GIVEN  BY  THE 
POINT  OF  TANGENCY  OF  THE  BOUNDARY  OF  THE  PRODUCTICN  FP^ONTIER 
AND  THE  HIGHEST  ATTAINABLE  INDIFFERENCE  CURVE,  WHERE  THE 

-h 

VECTOR  GF  THE  CONSUMER  IS  C  AND  THE  I-O  VECTOR  OF  THE  FIRM 

_f 

IS  Y  • 

BY  CONVEXITY  ASSUMPTIONS  THERE  EXISTS  A  SEPARATING  HYPER¬ 
PLANE  FOR  WHICH  THE  PRODUCTION  POSSIBILITIES  SET  LIES  ON 
ONE  SIDE  AND  THE  PREFERENCE  SET  ASSOCIATED  WITH  THE  HIGHEST 
ATTAINABLE  WELFARE  CURVE  LIES  ON  THE  OTHER  SIDE  OF  THE 
HYPERPLANE.  THE  HYPERPLANE  OR  PRICE  LINE  FUNCTION  : 


P  .Z  =  V 
WH  ER  E  J 

_•  — * 

V  =  P  .C  =  P  .Y 

IN  OUR  STUDY  WE  ASSUMED  THE  X  VECTOR  COMMODITIES  REFERS  TO 
GOODS  AND  SERVICES  SENSITIVE  TO  THt  OIL  POLLUTION  SUCH  AS 
FISHING  INDUSTRIES,  RECREATION  ACTIVITIES  ETC. 

ASSUMING  THAT  THE  PRODUCTION  POSSIBILITIES  SET  IS  A  FUNCTION 
OF  SPILL  SIZE.  WE  THEN  ASSUME  INCREASES  IN  POLLUTION  SIZE 
WOULD  CECPEASc  THE  PRODUCTION  OF  COMMODITIES  IN  THE  mREA 
FOULED  BY  OIL.  THIS  MAKES  SENSE  SINCE  FOULED  BEACHES  MEAN 
DECREASES  IN  INPUT  GF  RECREATION  PRODUCTION,  AND  FOULED 
WATERS  RESULT  IN  REDUCED  COMMERCIAL  FISHING  REVENUES. 

IF  GOODS  AND  SERVICES  AFFECTED  BY  THE  OIL  SLICK  HAVE  NO 
MARKET  VALUES  OR  PRICES,  W5  MAY  USE  A  SHADOW  PRICING  SCHEME 
TO  DETERMINE  THE  VALUES. 

THE  SHADOW  PRICE  OF  A  GIVEN  (CONSTRAINED)  COMMODITY  IS  DEFI¬ 
NED  AS  MARGINAL  VALUATION  OF  THE  COMMODITY  BY  RELAXING  THE 
CONSTRAINT  MEASURED  IN  TERMS  OF  THE  OBJECTIVE  FUNCTION  . 

THE  SOCIAL  COST  C(S)  OF  THE  OIL  POLLUTION  COULD  THUS  BE 
DESCRIBED  AS  THE  DIFFERENCE  IN  VALUE  BETWEEN  THE  PRODUCTION 
LEVEL  IN  THE  ABSENCE  OF  SPILLS  AND  THE  PRODUCTICN  LEVEL 
RESULTING  FROM  A  SPILL. 


C(S)  =  P(S=0) 


xl5=0)  -  X*S  =  1)  ) 
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VI.  COST  MODEL  DEVELOPMENT 


ANY  SPILLAGE  OF  OIL  INTO  THE  WATERS  REPRESENTS  A  LOSS  OF 
PROCLCTS  TO  THE  ECONOMY.  THIS  LOSS  IS  EQUAL  TO  THE  MARKET 
VALUE  CF  THE  PRODUCTS  LOST, 

IN  ADDITION  TO  THE  LOST  PRODUCT,  THERE  WILL  ALWAYS  BE  DAMAGE 
DONE  TO  THE  AQUATIC  ENVIRONMENT  AND  POSSIBLY  TG  THE  SURROUN¬ 
DING  LAND  AREAS  OR  BEACHES. 

MANY  FRACTIONS  OF  OIL  ARE  SOLUBLE  OR  EMULSIFIED  IN  THE  WATER 
AND  SINK.  SUCH  DISSOLVED  SUBSTANCES  CANNOT  BE  EASILY  TAKEN 
OUT  EY  SIMPLE  REMOVAL  OF  THE  VISIBLE  PRODUCT  AND  DAMAGES 
INFLICTED  MAY  NOT  BE  REVERSIBLE. 

THE  ENVIRONMENTAL  DAMAGE  FACTOR  IS  COMPLICATED  BY  THE 
EXISTENCE  OF  SHORT  AND  LONGER  TIME  EFFECTS. 

THE  LOCATION  OF  THE  SPILL  IS  CRITICAL  TO  THE  COST  THAT  IS 
INCURRED. 


A.  GENERAL  COST  MODEL. 


AS  WAS  PREVIOUSLY  MENTIONED,  THE  COST  OF  A  PARTICULAR  GIL 
SPILL  C(S)  : 


C  (  S  )  =  A(  S  ) 


WHERE: 

A  (S)  =  DIRECT  COST  OF  OIL  SPILL. 

B  (S)  =  INDIRECT  COST  OF  OIL  SPILL. 

8.  DIRECT  COST, 

DIRECT  COST  (S)  IS  EQUAL  TO  THE  MARKET  VALUE  OF  THE 
PRODUCT  LOSS 

A  (S)  =  C  .V 


V  =  VOLUME  OF  GIL  SPILLED 
C  =  MARKET  PRICE  PER  UNIT  VOLUME 
S  =  SPILL  INDICATOR. 


C.  INDIRECT  COST. 
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INDIRECT  COST  B(S),  IS  THE  VALUATION  OF  COMMOCITY  LOSS 
RESULTING  FROM  THE  OIL  SPILL. 


N 


B  (  S)  =  C.  Ax  =  ^  C;  . 


TO  SIMPLIFY  THE  COST  MODEL  TO  A  RESOLVABLE  LEVEL,  CONSIDER 
ONLY  COMMODITIES  THAT  MIGHT  HAVE  SIGNIFICANT  EFFECT  IN  THE 
CALCULATION  OF  COST. 

ASSUME  THE  SIGNIFICANT  INDIRECT  COST  WAS  THE  SUMMATION  OF 
LOSS  IN  WILD  BIRD  POPULATION,  LOSS  IN  FISH  POPULATION,  LOSS 
IN  PROPERTY  VALUES,  LOSS  TO  THE  BUSINESS  ACTIVITIES  IN 
COMMERCIAL  FISHING  AND  RELATED  INDUSTRIES  AND  LOSS  TO  THE 
RECREATION  ACTIVITIES  (  SPORT-FISHING,  BOATING,  TOURISM  ETC) 


THUS  IN  CUR  MODEL  COST,  N  EQUAL  TO  FIVE,  WE  DENOTE  : 


=  RELATIVE  PRICE  (WEIGHTING  FACTOR)  CF  A  BIRO. 
AXi  =  NUMBER  OF  BIROS  KILLED. 


THE  RELATIVE  PRICE  REFLECTS  THE  SHADOW  PRICE  WHICH  IS  PRE¬ 
SUMABLY  DETERMINED  UNDER  SOME  APPROPRIATE  CONSTRAINED  MAXIMI 
ZATICN  PROBLEM. 


C2  =  RELATIVE  PRICE  OF  A  FISH 
AX2  =  NUMBER  OF  FISH  KILLED 


AS  IN  C-1  ,USE  SHADOW  PRICE  TO  DETERMINE  THE  RELATIVE  PRICE 
C-.2  OF  A  FISH.  IN  THIS  CASS  THE  FISH  ARE  NOT  VALUED  AT 
THEIR  MARKET  PRICE  TO  AVOID  DOUBLE  COUNTING,  SINCE  THE 
EFFECT  ON  THE  MARKET  IS  COMPUTED  BELOW  AS  THE  REDUCTION  OF 
OUTPUT  IN  the  COMMERCIAL  FISHING  ACTIVITY.  THE  SHADOW  PRICE 
REPRESENTS  THE  LONG  TERM  EFFECTS  AND  ALSO  REPRESENTING  THE 
real  obvious  loss  in  the  NATURAL  RESOURCE  BASE. 


=  DROP  IN  MARKET  PRICE  OF  PROPERTY  VALUE  PER 
^  UNIT  AREA  RELATIVE  TO  THE  BASE  PRICE. 

AX5  =  AREA  OF  VALUABLE  B EAC HES /COAST  FOULED  BY  OIL. 

C-4  IS  A  FUNCTION  OF  LOCATION.  IN  OUR  STUDY  WE  DIVIDE  SAN 
FRANCISCO  BAY  AREA  INTO  FIVE  DIFFERENT  LOCATIONS.  THIS  CCULD 
BE  REFINED  INTO  ON  THE  SPOT  LOCATION  AS  PREDICTED  BY  THE 
SIMULATION  TECHNIQUE.  IN  REALITY  WE  SHOULD  ALWAYS  USE  MUCH 
FINER  GRID  SYSTEM  TO  OBTAIN  A  MORE  ACCURATE  COST  PICTURE. 


C4  =  MARKET  PRICE  PER  UNIT  PRODUCT  OF  COMMERCIAL 
FISHING  ACTIVITIES  (AND  RELATED  INDUSTRIES) 
RELATIVE  TO  THE  BASE  PRICE. 

/\X4  =  REDUCTION  ON  PRODUCTION  CAPACITY  OF  COMMER¬ 
CIAL  FISHING  ACTIVITIES  (AND  RELATED 
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INDUSTRIESJ 


C 


5 


MARKET  PRICE  PER  UNIT  PRODUCT  OF  RECREATION 
ACTIVITY  RELATIVE  TO  THE  BASE  PRICE. 

REDUCTION  IN  PRODUCTION  CAPACITY  OF  RECREA¬ 
TION  ACTIVITY. 


•D.  DETERMINISTIC  MODEL  TO  ESTIMATE  BIRD  OR  FISH  DEATH. 


TO  ESTIMATE  THE  BIRO  OR  FISH  KILLED,  WE  USE  DIFFERENTIAL 
EQUATION  AS  FOLLOWS: 


dX 

dt 


=.  -k.Y 


WHERE  : 

k  =  DEATH  RATE  OF  BIROS  OR  FISH  AS  A  FUNCTION  OF 
LOCATION,  SIZE  AND  TYPE  OF  CIL  SPILL. 

TYPE  OF  OIL  RELATES  WITH  ITS  TOXICITY. 

Y  =  AREA  SPILLED  AS  A  FUNCTION  OF  TIME. 

SOLVE  THE  EQUATION  BY  USING  EULER  APPROXIMATION  TECHNIQUE: 

NUMBER  OF  BIRDS  OR  FISH  KILLED  =  AX 

M 

Axl  5  k.At.  S  Y. 

lt=T  ^ 


E.  DETERMINISTIC  MODEL  TO  ESTIMATE  PROPERTY  DAMAGE. 


IF  THE  OIL  SLICK  LANDED  ON  THE  SHORE,  THUS  THE  SOCIAL 
COST  OF  DAMAGE  WOULD  BE  CONSIDERED  TO  BE  EQUAL  TC  THE  LOSS 
OF  VALUE  IN  MARKET  PRICE  OF  THE  PROPERTY  FOULED  BY  GIL. 


2 

StTh^-. 

2 

.A  Xj  ■s  C-  .TT. 


2 

TTPf  REPRESENT  THE  AREA  FOULED  BY  OIL.  Rf  MAY  BE  OBTAINED 
FROM  THE  SIMULATION. 
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F.  MODEL  TO  ESTIMATE  THE  SOCIAL  COST  OF  BUSINESS  ACTIVITY 


THE  STANDARD  PRODUCTION  FUNCTION  MODEL  IS  A  SUFFICIENTLY 
WELL  KNOWN  TOOL  OF  ECONOMIC  ANALYSIS. 

TO  MAKE  THE  MODEL  TRACTABLE,  WE  SIMPLIFY  BY  ASSUMING  THAT 
ONLY  A  SINGLE  FIRM  POSSESSES  ALL  BUSINESS  ACTIVITIES  IN  TH 
BAY  AREA, AND  THE  EXISTENCE  OF  SOCIAL  PREFERENCE  OVER  RECRE 
TION  ANC  COMMERCIAL  FISHING  ACTIVITIES. 

ASSUME  THE  PRODUCTION  POSSIBILITIES  SET  OF  THE  FIRM  IS 
KNOWN  AS  FOLLOWS  : 


WHERE  : 

X  =  PRODUCT  OF  RECREATION  ACTIVITY 

Y  =  PRODUCT  OF  COMMERCIAL  FISHING  (AND  RELATED  IN¬ 
DUSTRIES) 

a  ^  0  ,  b  ^0 

a  ,b  IS  A  FUNCTION  OF  SPILL  SIZE  AND  LOCATION. 


c  2 

r  ^  ■  f 

ao(A2^e  -H  6  ^  } 

-h  'F 

b  =  bo(^e 


(^1  =  THE  PROPORTION  OF  PRODUCTION  RELATED  WITH  THE 

^  BAY  ACTIVITIES  (  SWIMMING,  SPORT  FISHING  ETC) 

^2  =  DEGRADATION  CONSTANT 

^2  =  THE  PROPORTION  OF  PRODUCTION  NOT  RELATED  TO  THE 

^  BAY  ACTIVITIES  FOR  EXAMPLE  PEOPLE  IN  TRANSIT  ETC 

=  THE  PROPORTION  OF  PRODUCTION  WITH  LOCAL  INPUT 
RESOURCES. 

=  DEGRADATION  CONSTANT 

fX  =  THE  PROPORTION  OF  PRODUCTION  WITH  INTERLOCAL 
5  INPUT  RESOURCES. 

V  =  VOLUME,  SPILL  SIZE. 

eu  ,  =  MAXIMUM  PRODUCTION  LEVEL  FOR  EACH  ACTIVITY 

°  GIVEN  ALL  AVAILABLE  RESOURCES. 

a^,  bp,  5,6  IS  A  FUNCTION  OF  LOCATION. 
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THE  SOCIAL  PREFERENCE  FUNCTION  IS  GIVEN  AS 


W(s;  =  c.x'^.Y^ 


WHERE  : 

C  =  CONSTANT,  GIVEN  AS  A  FUNCTION  OF  LGCATIGN. 

C<,  :  KNOWN  AS  GIVEN  AN0CK+P<1»<^>0,  (3>0. 

BY  ASSUMING  THAT  THE  FIRM  TREATS  PRICES  AS  GIVEN  ANC  SEEKS 
TO  MAXIMIZE  ITS  PROFIT, 

MAXIMIZE  ITS  PROFIT, 


MAX  =  P  .X  +  P^  .Y 
X  A' 


SUBJECT  TO 


0 


BY  USING  LAGRANGE  MULTIPLIER  TECHNIQUE,  THEN 


PROFIT  IS  MAXIMIZED  IF  ; 


\ 


ON  THE  CCNSUMER  SIDE,  GIVEN  BUDGET  CONSTRAINT,  MAXIMIZE  HIS 
PREFERENCES 

MAX  W(S)  =  C.X  .Y 

SUBJECT  TO  Z  =  P  .X  +  P  .Y 

USING  LAGRANGE  MULTIPLIER  TECHNIQUE  : 


33 


PREFERENCeS  ARE  MAXIMIZED  IF 


o<  -z. 

_ 

EQUILIBRIUM  CONDITIQN  IF:  Z  =71=711 

THE  CONDITION  APPLIES  IF  : 


THE  SUPPLY  ANC  DEMAND  CURVES  OF  'X'  PRODUCT  ARE 


a'-p* 

^  TT 

V  c<rr* 

■ 


THE  SUPPLY  AND  DEMAND  CURVES  OF  'Y*  ARE  : 
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b^P 


* 

y 


TT  * 


IS.TT^ 
P*(^+  p) 


WHERE  : 


tt:*=  ( 


2*2 
a  + 


v2„*2 


1/2 


THE  PRICE  FUNCTION  IS  : 


=  P 


Y  .  a 
e 


2  2  1/2 
b(  b  -  ) 


G.  SENSITIVITY  ANALYSIS. 

AS  AN  EXAMPLE  WE  ANALYSE  THE  CHANGES  IN  SUPPLY  CURVES 
AS  WE  CHANGE  THE  PRODUCTION  POSSIBILITY  PARAMETER  : 

CHANGES  IN  SUPPLY  CURVES  : 


■aa 

Tia 


a.b^  (P  /P 
X  y 


a 


a 


a 


( 


2  2  2 
(P^/Py)^  +  b^  ) 

^  +  2b^  (Py/Px  ) 

2  2  2 
a  +  b  (Py/Px  ) 


3/2 

2 


4.0 


3  5 


CONCLUSICN 


1.  IF  'A*  DECREASES,  THE  SUPPLY  CURVE  OF  WILL 

SHIFT  TO  THE  LEFT,  IF  *A'  INCREASES  CONVERSE  HOLDS 

2.  IF  ‘A*  DECREASES  THE  SUPPLY  CURVE  OF  'Y*  WILL 
SHIFT  TO  THE  RIGHT,  IF  'A*  INCREASES  CONVERSE 
HOLDS. 


CHANGES  IN  DEMAND  CURVES  : 


SYo 


3a 


3  Vo  z.a.c<.P, 

j-  a  ' 


r^+  (3) 


CONCLUSION  : 

1.  IF  'A*  DECREASES,  BOTH  DEMAND  FOR  'X*  AND  'Y*  WILL 
SHIFT  TO  THE  LEFT. 

2.  IF  'A*  INCREASES,  CONVERSE  HOLDS. 

CHANGE  IN  THE  EQUILIBRIUM  POINT  : 
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=.  o 


^  Xc.  _  2.  <3  oC  ^  3^0^  \  '■^2. 

'2>a  ■  o<  + p,  ^  c?<4- 
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CQNCLUSICN 


1.  NO  CHANGE  IN  'Y'  PRODUCTION  LEVEL. 

2.  IF  'A'  DECREASES,  THE  PRODUCTION  LEVEL  OF  ‘X*  WILL 
DECREASE,  IF  ’A*  INCREASES  ,  CONVERSE  HCLDS. 

THE  SOCIAL  COST  OF  A  SPILL  IS  EQUAL  TO  THE  MARKET  VALUE  OF 
THE  FRGCUCT  LOSSES. 


SOCIAL  COST  =  Px.{Xg-  X'g  )  +  P^  .  (  Y^  -  Y^'  ) 


WHERE  : 

Px  AND  Py  REPRESENT  THE  PRICE  PER  UNIT  PRODUCT  OF  X 
/IND  Y, 

X„  AND  Y«  REPRESENT  THE  OPTIMUM  PRODUCTION  LEVEL 
OF  X  AND  Y  BEFORE  THE  SPILL  OCCURS. 

x'  AND  y'  represent  THE  OPTIMUM  PRODUCTION  LEVEL  OF 
0  a 

X  AND  Y  AFTER  THE  SPILL  OCCURS. 


SEE  FIGURE  8A,B  ;  FIGURE  9A,  ;  FIGURE  10A,B  FCR  GRAPHICAL 

ILLUSTRATION. 
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X  UNIT  PRODUCED 


-  PROD.  ^033  CURUe 

-  hPIK  PROFIT  LINE 
*  UELFRRE  CuRUE 


\  \  *  - 1 

fig-8b 

a  * 

K  UNIT  PRODUCED 


K  UNIT  PRODUCED 


*  SGPCLy  CiJRUE 


K  UNIT  PRODUCED 


V  UNIT  PRODUCED 


y  UNil  PROCUClO 


VII 


DATA 


ANALYSIS 


ASSUME  THE  PARAMETEi^S  OF  THE  COST  MODEL  IN  THIS  STUDY  ARE 
SHOWN  IN  TA6LE-2 


REGION 

PROD. PARM 

BEACH  VA 

,LU£ 

BIRD  PARM 

FISH  PARM 

SOUTHERN 

=40,000 

H=4,000 

N 

Cl=  2 

N 

£2=  2  N 

BAY  /  EAST 

=10,000 

L=2,500 

N 

k:^=0:01 

k2  =0.12 

SOUTHERN 

aQ  =60,000 

H=6,000 

N 

Ci=  3 

N 

€2=  2  N 

BAY  /  WEST 

bg  =40,000 

L=3,500 

N 

k;^=0.01 

k2  =0.10 

CENTRAL 

a^  =100,000 

H=7,000 

N 

Cl=  5 

N 

C2  =  2  N 

BAY 

bg  =  50,000 

L=3,500 

N 

ki=0.05 

k2  =0.0  5 

NORTHERN 

=20,000 

H=4,000 

N 

Ci=  3 

N 

C2  =  4  N 

BAY  /  EAST 

bo  =30,000 

L=l,500 

N 

0 

• 

0 

It 

r— 1 

.M 

k2  =  0 . 1 5 

NORTHERN 

aQ  =100,000 

H=8,000 

N 

Ci=  10 

N 

(^  =  1C  N 

BAY  /  WEST 

bo  =100,000 

L=5,000 

N 

=0.10 

k2  =0.15 

*  N=NUMERAIRE  ;  BASE 

PRICE,  ONE 

LIT 

ER  CIL  = 

ONE 

NUMERAIRE 

=!«  6EACF  VALUE  :  N  PER  HECTARE. 

*  BIRD  (FISH)  VALUE  :  N  PER  UNIT  KILLED. 

P4RM  =  PARAMETER 

*  H  a  L  :  H  STANDS  FOR  THE  VALUE  OF  BEACH  OR  OTHER 

property  PRIOR  TO  THE  OIL  SPILLAGE  AND  L 
AFTER  THE  SPILL,  THUS  THE  MARKET  VALUE  OF 
THE  LOSS  (C-5  )  IS  GIVEN  BY  : 

C5  =  H  -  L. 


TABLE-2 
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THE  BECPEiTICN  S  COMHSRCIIL  FISHING  COST  MODEL  PASArST'SRS; 


i, 

=  O.S 

i,  =  0.75 

4 

II 

O 

\ 

=  25X10"'^ 

0 

o 

It 

£3  =  0.25 

0( 

=  0.45 

fh  =  0.33 

p. 

=  150.  N 

p  =  VARIES  AS  A  FUNCTION  CF 

V 

COG 

=  COCROINAT 

latitude 

longitude 

c,  COMPUTED  IN  METERS  FROM  THE  ORIGIN 
"  ;  37»  25'  N 

:  122^  35’  W 

USING  T^-eSE  fpIf3Sj^.6,g3TA!NED^F|OH 

INDUSTRIES)  ACTIVITIti  TAeULAi-D 


TIGN 

(AND 


VALUES 

related 


IN  TAELE-5« 


LOSS  IN  RECREATION 

ACT.  AND 

COMMERCIAL  FISHING 

ACT. 

VOLUME  (LITE 

RS) 

LOCATION 

lOtOOO 

100,000 

1 9  00  C  t  0 0 c 

SOUTHERN 

BAY/EAST 

6,671. 

595,255. 

5,099,228. 

SOUTHERN 

bay/wsst 

10  ,005. 

392,382. 

7,648,843. 

CENTRAL 

EAY 

16,67 1. 

1,488,135. 

12,746,06C. 

northern 

6AY/HEST 

16  ,673. 

1 t433t  133* 

12,748,066. 

northern 

EAY/EAST 

3,335. 

297,626. 

2,549,613. 

TABLE 

-  3 

ACC lOEN 
GF  THE 


TS 

ACCI- 

WILL 


TO  OCT  ATM  Twc  ^nCTAL  COST  OF  ANY  PGLLUTION 

IRsiis  B  irfj!I!='SS"?KL!S»!'l.n.!":iS!iri 

SIX  =xil-PLE^CF°COMpiTING  THE  SOCIAL  COST  WAS  OETAINEO  AS 
FQLLCWS  : 


GVE  VEN' 
ON  Or 


45 


SAMPLES  OF  COMPUTING  THE  SOCIAL  COST  USING  THE  COMPUTER 
SIMULATION  TECHNIQUE: 


1.  COC  : 

(20  200  *,50700  ) 

^.2 

• 

0 

II 

HRS 

VOL  : 

100,000.  LTRS 

=  1.9 

HRS 

sj  = 

398  METERS 

=  5.8 

HRS 

LQCATICN 

AREA  FLD 

BEACH  FLD 

DEAD 

BIRD 

C5A0  FISH 

CENTRAL  EAY  4234.1 

49.3 

212 

635 

SOCIAL  COSTS  =  1,764,765.  N 


2.  CCO  :  (20200;50700)  T,2.  =  0.9  HRS 

VOL  :  L,OOC,OCO.  LTRS  T2.3  =  S.9  HRS 

:  944  METERS  =  13.3  HRS 


LCCATICN 

AREA  FLD 

BEACH 

FLO 

DEAD 

BIRD 

DEAD 

FISH 

CENTRAL  EAY 

6940. 

230. 

347 

1041 

SOCIAL  COSTS  = 

14, 731, 

877. 

N 

3.  COO  : 

(29000;64500 ) 

=  0.2 

HRS 

VOL  : 

10, 

COO  LTRS 

Tz3 

=  0.4 

HRS 

Rj  : 

168 

METERS 

=  lt8 

HRS 

LOCATION 

AREA  FLD 

BEACH 

FLD 

DEAD 

BIRD 

DEAD 

FISH 

NORTHERN 

ViEST 

EAY 

29.4 

NONE 

2 

C 

SOCIAL  COSTS  = 

26 

,753. 

N 
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,  coo  :  (29C00;64500 ) 
VOL  '  ItOOOfOOO  LTRS 
Rj  :  944  METERS 


Tgj  =8.9  HRS 
=  18.3  HRS 


location 


AREA  FLO  BEACH  FLO  CEAO  BIRD  CEAD  FISH 


JORTFERN  BAY  2089.6 
WEST 

SOCIAL  COSTS  = 

5.  CCC  :  {36000;21000  ) 

VOL  :  IfOCOtOCO  LTRS 
Rj  :  944  METERS 


78.3  209  313 


13»988»186.  N 

T,^  =  0.9  HRS 

T23  =8.9  HRS 

Tj  =  18.3  HRS 


LOCATION  AREA  FLD  BEACH  FLO 


SOUTHERN  BAY  2»432.  NONE 

WEST 


CENTRAL  BAY  llf705.5  280. 


DEAD  BIRO 


24 


1,171 


CEAO  FISH 


243 


11,706 


SOCIAL  COSTS  =  22,406,971.  N 


6.  COO  :  (37500;1S500  ) 
VOL  :  100,000  LTRS 
Rj  :  598  METERS 


T  =  0.4  HRS 
12 

T^3  =  1.9  HRS 
=  5.8  HRS 


LOCATION 


AREA  FLO  BEACH  FLD  CEAD  BIRO  CEAD  FISH 


SOUTHERN  BAY  1,640. 
EAST 


NONc 


16 


197 
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SOUTHERN  BAY 


3,168 


50 


32 


380 


SOCIAL  COSTS  =  1,714,799. 


THE  AGGREGATE  LOSSES  SHOW  THE  SIGNIFICANT  DIFFERENCE  BET¬ 
WEEN  SHALL  ANC  LARGE  SPILLS. 

THE  MOST  SIGNIFICANT  ECONOMIC  LOSSES  ARE  SUFFERED  BY  THE 
BUSINESS  ACTIVITIES  IN  THIS  ILLUSTRATION. 

NO  ATTEMPT  HAS  BEEN  MADE  IN  THIS  STUDY  TO  FORMULATE  THE  LONG 
TERM  SOCIAL  COST  OF  OIL  SPILLS. 

THE  MAGNITUDE  OF  COMPUTED  VALUES  IN  THIS  ILLUSTRATION  NATUR- 
RALLY  DEPEND  UPON  THAT  OF  THE  PARAMETERS,  WHOSE  VALUES  ARE 
ARBITRARILY  CHOSEN,  BUT  ILLUSTRATIVE  OF  POSSIBLE  PEAL  WORLD 
SITUATIONS. 
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VIII,  CONCLUSICN 


ThS  PURPOSE  OF  THIS  THESIS  IS  TO  DETERMINE  THE  SOCIAL 
COSTS  OF  OIL  SPILLS. 

THE  REALIZATION  THAT  SUCH  COSTS  DEPEND  ON  THE  ORIGINAL  LOCA¬ 
TION  ANC  SIZE  OF  A  SPILL  AND  THE  EVENTUAL  AREA  AFFECTED 
NECESSITATES  DEVELOPMENT  OF  A  MODEL  WHICH  PREDICTS  HOW  A 
GIVEN  OIL  SPILL  WILL  SPREAD  UNDER  PLAUSIBLE  CIRCUMSTANCES. 
FOR  THIS  WE  HAVE  SYNTHESIZED  THE  WORKS  OF  FAY  (REF-6},  SMITH 
(REF-4},  SIVADIER  AND  MIKOLAJ  (REF-1)  ,  CONOMOS  (REF-5}  BY 

INCORPORATING  WIND  (RANDOMIZED  DIRECTION  AND  VELOCITY}, 
ESTUARINE  NON  TIDAL  DRIFT  AND  EVAPORATION  FACTOR. 

IN  ORDER  TO  DETERMINE  SOCIAL  COSTS  OF  AN  OIL  S'^ILL  WE  ALSO 
NEED  TO  IDENTIFY  THE  COST  COEFFICIENTS  OF  FACTORS  DAMAGED 
SUCH  AS  WILDLIFE,  RECREATIONAL  ANC  COMMERCIAL  ACTIVITIES, 
ETC.  EXACT  DETERMINATION  OF  THESE  COEFFICIENTS  IS  BEYOND  THE 
SCOPE  OF  THIS  PAPER.  THUS  IN  THE  QUANTITATIVE  ANALYSIS 
GIVEN  HERE  FOR  ILLUSTRATION,  THEY  ARE  GIVEN  PLAUSIBLE,  BUT 
ARTIFICIAL  VALUES. 
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ooooooooooooooo 


APPEND!  X 


A 


COMPUTER  SIMULATION  PROGRAM  OF 
THE  SPREAD  AND  MOVEMENT  OF  OIL  SLICK 
IN  SAN  FRANCISCO  BAY 


//HANieiA  JOB  (2696t 1446,RL44) t 'MUDJI ARDJO  SMC  1374* 

//  EXEC  F0RTCLGP»REGIQN.GO=180K 
//FORT.SYSIN  CD  * 

THESIS  THESIS  THESIS  THESIS  THESIS  THESIS  THESIS  THESIS  TH 
THE  SPREAD  AND  MOVEMENT  OF  OIL  SLICKS  ON  THE  WATERS. 

OSL  =  OIL  SLICK  LEEWAY 

WIND  =  WIND  SPEED  MEASURED  AT  10  METERS  Ei-EVATIGN  IN  KNOTS 

AN  ESTIMATE  OF  WIND  DRIVEN  SURFACE  CURRENT  MAY  BE  OBTAINED 
FROM  THE  FOLLOWING  EQUATICN  PROPOSED  BY  TH0RAD5  (1914)  : 

sew  =  SURFACE  CURRENT  DRIVEN  BY  THE  WIND 
PHI  =  geographical  LATITUDE 
sew  =  0.0361  *  SQRT(WIND) /SORT ( (S!N( phi ) ) 

IF(  WIND  .GT.  11.64  )  SC W=0 .01 26*WIN0/SQRT ( ( SIN ( PHI ) ) 


DIMENSION  TIMES(80)t  RADIUSIBO) 

OIMENSI ON  X(  80) ,  Y(  80 ) 

DIMENSION  CCAA(900)t  00AA{900) 

DIMENSION  0SLV(80)»  SCWV(30),  WINDV(80) 

DIMENSION  DELT(  80),HTT(  80) 

DIMENSION  CC(  SO,  30)  ,  DD(  30,  80) 

CIMEMSION  CCA(  6400),  QOA (  6400) 

EGLIVALENCE(CCA,CC), ( CBA,0D) 

DIMENSION  XMV(80),  YMV(80) 

DIMENSION  XP(200),  YP ( 200 ) , XO ( 200 ) , YQ ( 2 00 ) 

DIMENSION  DZCIOO),  HZ(IOO),  VXL(80) 
dimension  XQ(200 ) , YQ( 200) ,XR(200) ,YR(200 ) 

DIMENSION  XS (200 ) , YS (200) , XT ( 200 ) , YT ( 200 ) 

DIMENSION  XU(200) ,YU(200) ,XV(200) ,YV(200) 

DIMENSION  AX (800 ) ,AY( 300) 

96  FORMAT(  6F5.1  ) 

97  FQPMAT(  2F10.1  ) 

98  FGRMAT(  F10.3  ) 

99  FaRMAT(  513  ) 

100  FCRMAT(/,10X,'M=‘  ,I3,3X,'N='  ,I3,3X,'KK='  ,I3,3X,'NK  =  ', 
1I3,3X,'KX=‘,  13, 3X, ' XMV( 1 )  =  '  , F3 . 1 , 3 X, • YM V ( 1 ) = '  ,  F3 . 1  ) 

101  FQRMAT(/,10X, 'VOLUME  OF  OIL  SPILLED  IN  L  IT E RS  = ' , F 1 2 . 1 ) 

102  FORMAT ( lOX, • RADIUS  OF  SLICK,  M',3X,'AT  TIMS,  HRS', 

23X, 'THICK  OF  0 1 L , CM ' , 3X , ' TH I CK  OF  GIL  P ART! CL ES , CM ' , 
33X, 'VOL  IN  LTRS' ,3X, 'NO: ' ) 

103  FORMAT ( 15X,F8.2, IIX, F6. 1 , 10X,F7.3, 15X,  F7.3  ,12X,F10.1, 
45X,I3  ) 

104  FORMAT (//, 18F8. 2  ) 

105  FCPMAT(  18F8.2  ) 

106  FORMAT (/, lOX, ' MAXIMUM  RADIUS  OF  OIL  SLICK  IN  METERS=', 
5F8.2  ) 

107  FGFMAT(///) 

108  FGRMAT(  •  4= 

109  FCFMAT(/,10X, '  CRITICAL  THICKNESS  IN  CM:',F7.3  ) 

110  FORMAT( /, lOX, ' THICKNESS  OF  MOL.  DIFFUSION,  CM=',F7.3) 

111  FQRMAT(/,10X,'  TIME  REQUIRED  AS  INFINITY,  HRS:',F6.1  ) 

112  FaRMAT(  '!'  ) 

113  FORMAT ( lOX, 'GAMA  1= ' ,F5.1,3X,'GAMA2=',F5.1,3X,'GAMA3  =  ', 
7F5. 1,3X, 'GAMA4=' ,F5.1,3X,'GAMA5=',F5.1,3X,' GAMA6=' , 
8F5.1  ) 

114  FORMAT( /, lOX, 'T12  =  ',F6.1,3X,'  T23  =  '  , F6. 1  ) 
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(^OOOOOOO  OOO  O  OOOoO 


\ 


2  01  FORMAT { lOX , F9 . 3 , 5X , F 9.3 , 3X , F 9. 3 , 2X , F9 . 2 , 4X , F9 . 2 , 6X ,  I  3) 
20  2  FORMAT  (  13X  , »  WI  NO  '  t  3X,  •  QSL  •  ,  9X,  •  ,  1  OX  ,  »  XM V •  ,  iOX  t 

1 ‘YMV* ,8X, ‘NO:*  ) 

203  FORMAT (  /  ) 

300  FORMAT  (  16F5.1  ) 

301  FGRMAT(  16F8.2  ) 

302  F0PMAT(  16F5.2  ) 

303  FGRMAT(/, lOX, '  RADIUS  OF  OIL  =‘,F8.2  t*  NZ=  ',14  ) 

304  F0RMAT(/,10X, 'TOTAL  AREA  FOULED  IN  HECT A R E= * , F 14 . 1  ) 

READ 

GAMM4  =  ANGLE  BETWEEN  THE  WIND  VECTORS  AND  THE  COORDINATE 
SYSTEM 

READ(5,96)  GAMAlt  GAMA2t  GAMAS t  GAMA4,  GAMA5,  GAMA6 
SPILL  COORDINATE  MEASURED  FROM  ORIGIN  IN  METERS. 
RFAO(5,97)  XMV(L),  YMV(l) 

GIVE  THE  INPUT  DATA  OF  THE  QUANTITY  OF  OIL  SPILLED  IN  KILO 
LITERS. 

1  KILO  LITER  =  1000  LITERS. 

REAO(5,98)  QTT 

M,NtKK  =  FIXED  NUMBER  OF  COORDINATE  JOINTS  CF  SAM  FRAN 
CISCO  BAY  AREA  AND  CURRENT  VECTORS  ON  THE  BAY  AS  PRECICTED 
BY  ESTUARINE  NONTIDAL  DRIFT  TO  BE  PLOTTED  BY  THE  CGMOUTER. 
NK  ;  NUMBER  OF  ITERATION  RELATED  TO  TIME  NEEDED  TO 
ESTIMATE  THE  OIL  MOVEMENT. 

KX  :  NUMBER  OF  ITERATION  TO  COMPUTE  THE  COORDINATE  OF 
THE  OIL  SPREAD. 

LIMIT  OF  (NK*KX)  =  900  POINTS. 

REA0(5t99)  M,  N,  KK,  NK ,  KX 
C  COORDINATE  POINTS  TO  CRAW  SF  BAY  AREA. 

R5AD(5t300)  (  XO (  I  )  ,  1  =  1, KK  ) 

READ(5,300)  (  YO(I  ),  1=1, KK  ) 

C  THE  COORDINATE  POINTS  TO  DRAW  THE  ESTUARINE  NON  TIDAL 
C  DRIFT. 

READ(5,302)  (AX(I),I=1,N  ) 

REAOt  5, 302)  (AY( I ) ,I  =  1,N  ) 

C  the  COORDINATE  POINTS  TO  DRAW  THE  ESTUARINE  NON  TIDAL 
C  DRIFT  CONTINUED 

R6AD(  5,  202)  ( XR ( I  )  , I  =  1 , M) 

REA0(5,202)  ( YR ( I  )  , I  =  1 , M ) 

0444444444444444444444444444444444444444444444444444444 444 44 
C 


C 

c 


Cl  =  0.53 
C2  =  0.026 
ZU  =  0.01 
RhCIL  =  0.95 
RHOWAT  =  1.0 
GRAV  =  980.0 
D  =  1.0 
SIGMA  =10. 

PHI  =  38. 

IX  =  1234567 

OELRHG=  RHOWAT  -  PHOIL 
OELTA=  DELRHO/RHOIL 
VOL  =  GTT=^  1000000.0 
VOLUME  =  VQL/1000. 


COMPUTE 


C  COMPUTE  THE  CRITICAL  THICKNESS  OF  THE  OIL  SLICK. 

HC  =  (  SIGMA/(OELRHO  #  GRAV))«-.5 
WRITE! 6, 112) 

WRITE(6, 108) 

WRITE(6,100)  M,  N,  KK,  NK ,  KX,  XMV(l),  YMV(l) 
WRITE(6,113)  GAMAl,  GAMA2,  GAMA3,  GAMA4,  GA^AS,  GAMA6 
WRIT£(6,101)  VOLUME 
WRITE(6,109)  HC 

C  MAXIMUM  RADIUS  OF  GIL  SLICK  AS  ONE  UNIT  AT  ZERO  GROWS  RATE 
C  (  INDEPENDENT  OT  TIME  ) 
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VVV  =  VQLUM=/1000. 

REND  =  100.  *  (  10.=!'*5*VVV*=!=0.  75/  3.14  )^*0.5 
CALL  RANDU(IX,IY, YFL) 

C3  =  0.78  +  YFL*0.02 
C  ENC  RADIUS  REDUCED  DUE  TO  EVAPORATION 
REND  =  RENO  *  C3  /2. 

TINF=  (  REND/2.3  )=^*4*RHCWAT**2*ZU/S  IGMA^^Z 

TINF  =  TINF*=i'0. 33333 

HINF  =  SQRT(  TINF*D/10.*=S‘5 ) 

TINF  =  TINF/3600. 

REND  =  RENO/100. 

T12  =  (  1.  45/1.1 4  VOL /(DELTA’f'GRAV^ZU)  )  =!'*0.333 

TL2  =  T12/3600. 

T23=(  1.45/2.  30)*=«'2=!'(RH0WAT/SIGMA  1  *  VOL=;'*0. 6667 
1(  DELTA  *  GRAV  ZU  )**0.333 
T23  =  T23/3600. 

VVL  =  VOL 
WRITE(6,106)  REND 
WRITE(6,110)  HINF 
WRITE{6tlll)  TINF 
WRITE(6,114)  T12,T23 
WRITE(6t 108) 

WRITE! 6,107) 

WRITE(6,102) 

TIMES! 1 )  ^  0. 

TS  =  0. 

DC  20  K  =  1,NK 

C  ASSUME  THE  WIND  BLOWS  RANDOMLY  AT  SPEED  BETWEEN  0.5  TC  25 
C  KNCTS  PER  HOUR. 

CALL  RANDU!! X, lY, YFL) 

WIND  =  YFL  24.5  +  0.5 
XXX  =  SIN!0.01745*PHI ) 

C  COMPUTE  THE  SEA  SURFACE  CURRENT  VELOCITY  INDUCED  BY  THE 
C  WINDS. 

sew  =  0.0361  *  SQRT!WIND)/SQRT!XXX) 

IF!  WIND  .GT.  11.64  )  SCW  =  0 . 0 1 26*WI ND/ SQR T ! XXX ) 

C  COMPUTE  THE  OIL  SLICK  LEEWAY  BLOWS  BY  THE  WIND. 

OSL  =  0.0199  =>=  WIND 

IF!WIND.GT.5.0.AND.WIND.LT .25.0 )  OSL= . 0 1 79* wl NO+0 . 0 196 
C  CONVERT  KNOTS/HR  TO  METERS/HR  : 

WINDV!K)=  1850.  *  WIND 
SCWV!K)  =  1850. *SCW 
OSLV!K)  =  1350. *OSL 
C  SPREAD  BY  GRA V IT Y- INE RTI A  FORCE 
TIME  =  FLOAT!K) 

T  =  TIME*15,*60. 

TIME  =  TIME/4. 

IF  (TIME  .GE.  T12)  GO  TC  1 
DELT!K)  =  !ZU*T)**.5 
VLL  =  vaL/10.**4 

R1  =  10.*!  DELTA  *  GRAV  *  VLL  *  T**2  )**0.25*1.14 
HT  =  V0L/R1**2 
R1  =  Rl/100. 

RADIUS!K)  =  R1 
TS  =  TS  +  1. 

GO  TO  21 

C  SPREAD  BY  GRAVITY-VISCOUS  FORCE 
1  AA  =  ZU**.5 

TIME  =  FLOAT!K)  -  TS 

IF!TIME  .GT.  24.)  TIME=!TIME  -  24.)*2.  +  24. 

IF!TIM£  .GT.  48.)  TIME=!TIME  -  48.)*3.  +  48. 

T  =  TIMe*3600. 

IF  !TIME  .GE.  T23)  GO  TO  2 
0ELT!K)  =  !ZU*T)**.5 

R 2= 1.45*! ! !OELTA*GRAV*VaL**2*T**1.5)/AA ) **.166667 ) 

HT  =  V0L/R2**2 
R2  =  R2/100. 

RADIUS!K)  =  R2 

IF!  RADIUS!K)  .GE.  RENO  )  GO  TO  9 
GO  TO  21 

C  SPREAD  BY  SURFACE  TENS IGN-VI SCOUS 
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2  65=  ZU*RHCIL*^2 

DtLT(K)  =  { ZU’f'T)  *=i=.5 

R3  =  ((SIGMA^^PZ  =>=  T**3  )/  BB  \**.Z5  ^  2.20 

HT  =  V0L/R3=«'=!'2 
R2  =  R3/10Q. 

RACIUS{K)  =  R3 

IF(  RADIUS(K)  .G6.  REND  )  GO  TO  9 
GO  TO  21 

9  RflCIUS(K)  =  REND 

HT  =  vaL/(REND*100.)*^‘2 
C 

C  LOSS  CUE  TO  EVAPORATION 

21  FF  =  C1=!'TIME=«'60./ (  1.  +  C2’!‘TIME=!'60.  ) 

VOL  =  { VVL  -  VVL*FF*0.01) 

VXL(K)  =  VOL/1000. 

C 

C  PRINTS  THE  RADIUS  OF  OIL  SPREAD  IN  NETERS,  THE  TIME  OF 
C  SPREAD  IN  HOURS,  THE  THICKNESS  OF  OIL  AS  A  CUNCTICN  OF 
C  TINE  IN  CM,  THE  THICKNESS  OF  OIL  LAYER  IN  CM,  AND  THE  RE- 
C  MAINING  VOLUME  OF  OIL  AFTER  EVAPORATICN  IN  LITERS. 

HTT(K)  =  HT 

WRITE (6,103)  RAOIUS{K) ,TIME ,HTT(K) ,0£LT(K) ,VXL{K) ,K 
X(K)  =  ALOGIO(TIME) 

Y(K)  =  AL0G10(RADIUS(K) ) 

ALPHA  =  0.0 
TINES(K+1)  =  TIME 

C  THE  ESTUARINE  NCNTIOAL  DRIFT  WAS  SIMULATED  AS  A  UNIFORM 
C  CURRENT  WITH  CONSTANT  SPEED  APPROX.  65  M/HOURS  IN  THE 
C  SOUTH  BAY  AND  200  M/HOURS  IN  THE  NORTHERN  BAY  AREA. 

CRAN  =  62.500«(TIMES(K+1)  -  TIM£S(K)) 

IF(  XMV(I)  .GE.  16125.  .AND.  YMV(I)  .GE.  51600.  ) 
ICRAN  =  208.3*{TIM£S(K+1)  -  TIM£S(K)) 
l^ITIMES (K+l)  -  TIM£S(K)) 

BETAl  =  0.01745*( YFL^IOO.  +  90.  ) 

IF{  XMV(K)  .GT. 35475.  .AND.  YMV ( K  )  .GT . 806 3 . 0  ) 

13ETA1  =  0.01745*(  YFL^IOO.  +  65.  ) 

IF(  XMV(K) .GT. 21736.  .AND.  YMV (K )  .GT.  17737.  ) 

2BETA1  =  0.01745=«'(  YFLY60.  +  65.  ) 

IF(  XMVIK)  .GT.21  736.  .AND.  YMV ( K ) .GT . 3 37 00 .  ) 

38ETA1  =  0.01745^(  YFL=^100.  +  45.  ) 

IF(  XMV(K),GT.  16125.  .AND.  YM V ( K  )  .GT .5 1 600  .  ) 

4BETA1  =  0.01745*{  YFL't'60.  +  160.  ) 

IF(  XMV(K)  .LE. 16125.  .AND.  YMV ( K  )  . LE . 51 60C . 0  ) 

5BETA1  =  0.01745*(  YFL^130.  +  90.  ) 

XCRAN  =  CRANYCOS (BETAl ) 

YCRAN  =  CRANES  IN ( BETAl  ) 

IX  =  lY 

CALL  RANDUdX,  lY  ,YFL  ) 

BETA  =0 .01745*( YFL^GAMAl  +  GAMA2  ) 

IF(  XMV(K)  .LE.  16125.  )  6ETA=  0 . 0  1745=!=  (  YFL  ^1  80  .  +  9  0  .  ) 
IF(  XMV(K)  .GT.  16125.  .AND.  YMV(K)  .GT .  36230.  ) 
IBETA  =  0.01745*(  YFL=!=GAMA3  +  GAMA4  ) 

IF(  XMV(K)  .GT.  16125.  .AND.  YMV(K)  .GT.  51600.  ) 
2BETA  =  0.01745^(  YFL=i‘GAMA5  +  GAMA6  ) 

ROUTE  =  0SLV(K)*(TIMES(K+1)-TIM£S(K) ) 

XMV(K  +  1)  =  XMV(K)  +  ROUTE=^COS  (BETA)  +  XCRAN 
YMV(K  +  1)  =  YMV(K)  +  ROUTE=!=  S I N  (  B  ET  A  )  +  YCRAN 
IX  =  lY 

C  GENERATE  RADIUS  OF  OIL  SPREAD. 

DC  30  L  =  1,KX 

ALPHA  =  ALPHA  +  360 . / FLO AT ( KX ) 

CC(K,L)  =  XMV(K+1)  +  RADIUS  (K)=!=CGS(  ALPHA  ) 

CC(K,L)  =  YMV(K  +  L)  +  RADI  US ( K ) *S I N ( ALPHA  ) 

30  CONTINUE 

IF(  HTT(K)  .LT.  0.010  )  GO  TO  22 
20  CCNTINUE 

22  WRITE(6,112) 

WPITE(6,202) 

CO  50  K=1,NK 

WRITE(6 , 201)  WINOV(K) , OSLVi K) ,3CWV(K) , XMV( K+1) , 
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2Y^V( K+1 ) tK 
50  CONTINUE 

V^PITE(6,112) 

C  PICT  THE  SLICK  MOVEMENTS  INDUCED  BY  THE  WIND  ftNC  ITS 
C  SPREAC. 

NN  =  NK^KX 
L  =  0 

DC  999  11=1, NK 
OG  999  JJ=1,KX 
L  =  L  +  L 

C  LET  CONVERT  50000  METERS  EQUAL  TO  10  INCHES: 

CCA(L)  =  CC(  II, JJ  ) 

CCAA(L)  =  0.00620*(CCA(L)  ) 

IF(  CCAA(L)  .LE.  O.O  )  CCAA(L)  =  O.O 

IF(  CCAA(L)  .GE.  9.0  )  CCAA(L)  =  9.0 

CCA(L)  =  DD( I  I, JJ  ) 

DCAA(L)  =  0.00020=!'DDA(L  ) 

IF(  DDAA(L)  .LE.  0.0  )  OOAA(L)  =  0.0 

IF(  CCAA(L)  .GE.14.0  )  DOAA(L)  =  14.0 

999  CONTINUE 
C 

C  PLOT 

C 

C  PLCT  IN  A  LIN-LIN  SCALE  THE  THICKNESS  CF  OIL  AS  A  FUNCTION 
C  OF  TIME 

CALL  PL0TP(X,HTT,NK,0) 

WRITE( 6, 112) 

C  PLCT  IN  A  LIN-LIN  SCALE  VOLUME  AS  A  FUNCTION  CF  TIME. 

CALL  PLCTP(  X,VXL,NK,0  ) 

V<RITE(6,1121 

C  PLCT  IN  A  LOG-LOG  SCALE  THE  RADIUS  OF  SPILL  AS  A  FUNCTION 
C  OF  TIME 

CALL  FL0TP(X,Y,NK,0) 

C 

REAL  XTITLE( 2 )/' LATI ' , 'TUDE '/ 

REAL  YTITLE(  3)  /  '  LONG'  ,  •  ITUD' ,  '  £  ' /’ 

XI  =  0.0 
X2  =  0.0 
NCX  =  -8 
SIZEX  =  8. 

XMIN  =  0.0 
CX  =  5C00. 

Y1  =  0.0 
Y2  =  0.0 
NCY  =  9 
SIZEY  =  14. 

YMIN  =  0.0 
OY  =  5000. 

DRAW 


C 


NK  =  NK+1 
CALL  PLCTS 

CALL  AXIS( XI, X2, XTITLE, NCX, SIZEX,  0.0, XMIN 
CALL  AXIS( Y1,Y2, YTITLE,NCY, SIZEY, 90. , YMIN, 
NZ  =  0 

CO  911  1=1, NK 

HTT(I)  =  O.QOOZO^i  XMV(I)  ) 

C£LT(I)=  0.0002-YMV(I ) 

NZ  =  NZ+1 

IF(  HTT(I)  .LT.  0.0  .OR.  OELT(I)  .LT.  0.0 
IF(  HTT ( I ) .GE.8.0  .OR.  DELT ( I ) . GE  .  14.0 )  GG 
911  CONTINUE 

7  CALL  LIN£(HTT,0£LT,NZ,  1,  2  ) 

CALL  LINE (CCAA,DDAA,NN,  1,-7) 

CALL  PLCT  (0.0, 16. 0,-3 ) 


t  DX) 
CY  ) 


)  GO  TO  7 
TO  7 


IF(  YMV(l)  .GE.  36757.5  )  GO  TO  555 
IF(  XMV(l)  .GE.  48165.  )  GO  tq  556 
IF(  YMV(l)  .GE.  60850.  )  GO  TO  557 


55 


ooo 


IF(  XMVd  )  .GT. 58000.  .OR.  YMV  (  1  )  .GT.  7C900.  J  STOP 

SOUTHERN  BAY 

DC  912  I=ltKK 
XC(I )  =  X0(I )  -  140. 

!F(  xg(  I  )  .L£.  0.0  )  XQd  )  =  0.0 
XG(I)  =  XQd  )4:0. 0725 
YQ(I  )  =  YOd  ) 

IF(  YQ(I)  .GE.  125.  )  YQd  )  =  125. 

YQd)  =  YQd  )=i'0. 0725 

912  CONTINUE 
SUM  =  0.0 
NZ  =  0 

CO  401  I=1,MK 

HTTd)  =  0.00039^XMV(I  )  -  14.0 
CELTd)  =  0.00039*YMV(  I  ) 

IF(HTT(I).LT.0.0.  OR . DELT ( I ) . LT . 0 .0 )  GO  TO  401 
IF  (HTT(  I  )  .GT  .9.0  .OR.  OELTd  )  .GT.  12.0)  GO  TO  401 
CZd)  =  (  YMV(I)  -  YMV(I-l)  )^^2 
HZd)  =  (  XMVd)  -  XMV(  I-l  )  )*=!=2 
NZ  =  NZ+1 

FA  =  2.*RA0IUS(I  j^^'SQRTi  DZd)  +  HZ(I)  ) 

SUM  =  SUM  +  FA 

IF(  CELTd)  .LT.  3.5  )  GO  TO  666 
C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

HTTd)  =  hTT(  1)^0.725 
CELTd)  =  OELTd  )*0.725 

IF(  HTT  (  I  )  .GE.3.0  .OR.  CELTd). GE.  8.4  )  GO  TO  3 
401  CONTINUE 

8  IF(  NZ  .LT.  3  )  GO  TO  888 

CALL  LIN£(  HTT, CELT, NZ,  1,  2  ) 

NM  =  NZ+KX 
CC  501  I=1,NM 

CCAAd)  =  0.00029=i'CCA(  I )  -  14. 

IF(  CCAA(I)  .LE.  0.0  )  CCAAd)  =  0.0 
IF(  CCAAd)  .GT.  4.0  )  CCAAd)  =  4.0 
OCAA(  I)  =  0.00039’!‘DOA(  I ) 

IF(  DOAA(I)  .GE.12.0  )  ODAAd)  =  12.0 
C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

CCAAd  )  =  CCAA  ( I  )  *0.725 
.  CCAA(  I)  =  CCAAd  )*0.725 
501  CONTINUE 

SLMl  =  0.0001*SUM 
WRITE(6,304)  SUMl 
WRITE(6,303)  RAOIUS(NZ),  NZ 
CALL  LINS(  XQ, YQ,KK, 1,1) 

CALL  SYMBOL(  2.3,  9. 4, . 14, » EAST  SOUTHERN  BAY', .0,17  ) 
CALL  LINE(CCAA,OOAA,NM,  1,-7  ) 

CALL  PLOT  (0.0, 12. 0,-3 ) 

C 

838  DC  913  1=1, KX 

XS(  I  )  =  XOd  )  -  60. 

IF  (  XS( I )  .LE.  0.0  )  XS( I )  =  O.C 
IF  (  XSd  )  .GE.  90.  )  XSd  )  =  90. 

XSd  )  =  XSd  )*0.0725 
YSd  )  =  YOd  )  -  50. 

IF(  YSd  )  .LE.  0.0  )  YSd  )  =  0.0 
IF(  YSd)  .GE.  120.  )  YSd)  =  120. 

YSd)  =  YSd  )*0.0725 

913  CONTINUE 
SUM  =  0.0 

HTTd)  =  0.00039*XMVd)  -  6.0 
DELT(l)  =  0.00039*YMV(  1)  -  5.0 

IF(  HTTd).Gc.  0.0  .AND.  OELTd). GE.  0.0  )  NZ  =  0 
CO  701  1=1, NK 

HTTd)  =  0.00039*XMVd  )  -  6.0 

OELTd)  =  0.00039*YMVd  )  -  5.0 

IF (HTT (I  ).LT.O.O.  OR . CELT (I ) . LT . 0 . 0 )  GO  TC  701 

IF(HTT( I  ).GT.9.0.  OR. CELT (I ) .GT . 12.0 )  GO  TO  701 

CZd)  =  (  YMV(I)  -  YMVd-l)  )**2 

HZd)  =  (  XMVd)  -  XMVd-l))**2 


56 


oooo 


NZ  =  NZ+l 

F4  =  2.*RA0IUS(I  )*SQRT(  DZ ( ! )  +  HZ(I)  > 

SUM  =  SUM  +  FA 

IF(  HTT(I  )  .GT.5.25.AiN0.DELT(I  )  .GT.9.0  )  DEL7(I) 

IF(  HTT(I),GT,8.0.AN0.DELT(I ).GT.6.5  )  OELT(I)= 

C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

FTT(I)  =  HTT(I)«0.725 
DELT( I )  =  D£LT( I )*0.725 
IF{  HTT(I)  .LT.  2.4  .ANO.DELTd)  .LT.  7.3  )  GO  TO  3 
IF(HTT( I).GE,5.5  .AND.  OELT ( I ) .GE . 5. 0  )  GO  TO  3 
701  CONTINUE 

3  CALL  LINF(HTT,DELTtNZT  It  2  ) 

NM  =  NZ^KX 
CO  601  I=1,NM 

CCAA(I)  =  0.0C03  9=!'CCA(  I  )  -  6. 

IF(  CCAA(I)  .L£.  0.0  )  CCAA(I)  =  0.0 
IF(  CCAAd)  .GE.  9.0  )  CCAAd)  =  9.0 
CCAA(I)  =  0.00039*D0Ad  )  -  5.0 
IF(  CCAAd)  .GE.  12.0  )  OOAA(  I  )  =  12.0 
IF(  CCAA(  n.GT.5.25.AN0.D0AAd  )  .GT.  9.0)  DDAAd)=  9.C 
IF(  CCAAd  ).GT.3.0.ANO.DOAAd)  .GT.6.5  )  OOAAd)=  6.5 
C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

CCAAd  )  =  CCAAd  )*0.725 
ODAA(  I )  =  DOAAd  )  *0.725 
601  CONTINUE 

SLM2  =  0.0001*SUM 
WRITE(6»304)  SUM2 
WPITE(6t303)  RADIUS(iMZ),  NZ 

CALL  SYMBOL(  2.3,  9. 2 , . 14, • WEST  SOUTHERN  BAY', .0,17  ) 
CALL  LINE(  XS,  YS,KK,  1,  1  ) 

CALL  LINE(CCAA,DOAA,NM,  1,-7  ) 

CALL  PLOT  (0.0, 12. 0,-3) 


NORTHERN  BAY 

556  DC  916  1=1, KK 

XVd)  =  XO(I)  -  140. 

IF(  XV(  I  )  .LE.  9.0  )  XVd  )  =  0.0 
XVd)  =  0.0725*XVd) 

YVd  )  =  YOd  )  -  230. 

IF(  YV(  I  )  .LE.  0.0  )  YVd  )  =0.0 
YVd)  =  0.0725*YVd) 

916  CONTINUE 
SUM  =  0.0 
NZ  =  0 

DC  705  1=1, NK 

HTTd)  =  0.00039*XMV(  n  -  14.0 
OELT(I)  =  0.  00039*YMVd  )  -23,0 
IF (HTT( I ) .LT .0.0 .  GR . OELT (I ) .L T . 0 .0 )  GO  TO  705 
DZd)  =  (  YMVd)  -  YMVd-1)  )**2 
HZd)  =  (  XMVd)  -  XMV(I-1))**2 
NZ  =  NZ+L 

FA  =  Z.TRADIUSd  )*SQRT(  DZd)  +  hZ(I)  ) 

SUM  =  SUM  +  FA 
C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

HTTd)  =  HTTd)*0.725 
DELT( I )  =  OELTCI )*0.725 

IF(  HTTd). GT.6.5  .OR.  CELT  d  ).  GT  .  5.0  I  GO  TO  6 
705  CONTINUE 

6  IF(  NZ  .LT.  3  )  GO  TO  666 

CALL  LIN£(  HTT,DELT,NZ,  1,  2  ) 

NM  =  NZ*KX 
DO  602  1=1, NM 

CCAAd)  =  0.00039*CCA(I  )  -  14. 

IF(  CCAAd)  .LE.  0.0  )  CCAAd)  =  0.0 

IF(  CCAAd)  .GE.  9.0  )  CCAAd)  =  9.0 

DCAA(I)  =  0.0C039*D0A d )-  23. 

IF(  DCAA(I)  .LE.  0.0  )  DDAAd)  =  0.0 

IF(  DDAA(I)  .GS.  7.0  )  DDAA(I)  =  7.0 

C  CONVERT  TO  THESIS  FORMAT  3.5X11. 

CCAAd)  =  CCAAd  )*0.725 
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CC^A(!)  =  DDAA(  I  )=!'0.725 
602  CONTINUE 

SUMS  =  0.0001*SUM 
WRITE(6,304)  SUMS 
WRITE{6,303)  RADIUS(NZ),  NZ 

CALL  SYMBOL(  2 . 1 , S .43 , . 14 , • EAS T  NORTHERN  eAV',.0, 
CALL  LINE<  XV,YV.KK,  1,  1  ) 

CALL  LINE(CCAA,DDAA,NM,  1,-7 


CALL  PLQT(  0.0, 10.0, -3  ) 


) 


XU(  I  ) 
XU(  I  ) 


0.0 

90. 


YU(I )  =0.0 


SS7  DO  91S  1=1, KK 

XU(  I  )  =  XOd  )  -  SO. 

IF (  XU( I  )  .LE.  0.0  ) 

IF(  XU(  I  )  .GE.  90.  ) 

XU(I)  =  0.0725-XU(I) 

YU(I)  =  Y0(!)  -  230. 

IF(  YU( I )  .LE.  0.0  ) 

YU(I)  =  0.0725-YU(I) 

91S  CONTINUE 
SUM  =  0.0 

HTT(l)  =  0.00039=f=XMV(  1)  -  S.O 
D£LT(L)  =  0.00039*YMV(  1)  -  23.0 
IF(  HTT(1).GE.  0.0  .AND.  DELT(l) 

DO  704  1=1, NK 

hTT(I)  =  0.00039=i'XMV(  I  )  -  S.O 
\  DELT(I)  =  0.00039*YMV( I )  -25.0 
IF(  HTT(I)  .LE.  1.0  )  HTT(I)=1.0 
IF(  OELT(I).G£.  S.8  )  0ELT(I)=S.8 
IF(HTT( I ) .LT.0.0 .  OR.DELT( I  ).LT  .0.0) 

OR. OELT ( ! ) .GT .  6.0)  GO 


GE.  0.0  )  NZ  =  0 


6.3,  .14,  'CENTRAL 

-  YMV(I-l)  )’i'’*'2 

-  XMV( I-l) 


)*3QRT(  DZd  )  +  hZ(  I)  ) 


IF(HTT( I ) .GT.9.0  . 

CALL  SYMBOL(  2.2, 

DZd)  =  (  YMV(I) 

HZd)  =  (  XMV(I) 

NZ  =  NZ+1 
FA  =  2.  ’^RADIUSC  I 
SUM  =  SUM  +  FA 
C  CONVERT  TO  THESIS  FORMAT  8.SX11. 
HTTd  )  =  HTT(  I  )*0.72S 
CELT( I )  =  DELT( I )*0.72S 
IF(  HTTd)  .LE.  0.72. OR.  OELT(I) 
704  CONTINUE 

S  CALL  LINE(  H7T,DELT,NZ,  1,  2  ) 

NM  =  NZ*KX 
DO  604  I=1,NM 

CCAA(I)  =  0.00039=«=CCA(  I  )  -  S. 

IF(  CCAA(I)  .LE.  1.0  )  CCAAd)  = 

IF(  CCAA(I)  .GE.  9.0  )  CCAAd)  = 

OCAA( I)  =  0.00039*DDA( I )-  23. 

IF(  DDAA(I)  .LE.  0.0  )  DDAA(I)  = 

IF(  OCAA(I)  .GE.  5.3  )  DDAA(I)  = 

C  CONVERT  TO  THESIS  FORMAT  8.5X11. 
CCAAd  )  =  CCAAd  )=^0.725 
CDAA(  I)  =  DDAAd  )’«'0.725 
604  CONTINUE 

SUM4  =  0.0001*SUM 
WRITE! 6,304)  SUM4 
WRITE(6,503)  RADIUS(NZ), 

CALL  LINE!  XU,  YU,KK,  1, 


BAY' 


TO  704 

0,  11  ) 


.GE.  3.75)  GO  TO 


l.C 

9.0 


0.0 

5.3 


NZ 
1  ) 


CALL  SYMBOL!  2 . 1 , 5 .43 , . 14 , • WEST  NORTHERN  BAY'r.O, 
CALL  LINE(CCaA,DDAA,NM,  1,-7  ) 

CALL  PLOT!  0.0, 10.0, -3  ) 

CENTRAL  BAY. 


555  DC  914  1=1, KK 

XT!I )  =  XO! I )  -  40. 
IF!  XT! I )  .LE.  0.0  ) 


XT!! )  =  0.0 


IF! 

XT!  I  ) 

.GE. 

90 .  ) 

XT! I )  =  90. 

XT!I)  =  0 

.0725*XT!I  ) 

YT!I )  =  YO!I  ) 

-  145. 

IF! 

YT!  I  ) 

.LE. 

0.0  ) 

YT! I )  =0.0 

IF! 

YT!  I  ) 

•  GE  • 

90.  ) 

YT!I )  =  90. 

17  ) 


5 


17  ) 
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YTU»  =  0.0725*Y7(I) 

914  CCNTINUE 
SUM  =  0.0 

HTT(l)  =  0.00039=!'XMV(  1)  -  4.0 
OeLT(l)  =  0.00039*YMV( 1)  -  14.5 

IF(  HTT(1).GE.  O.O  .ANO.  DELKD.GE.  0.0  )  NZ  =  0 
00  703  1=1, NK 

hTT(I)  =  0.00039*XMV( I  )  -  4.0 
OELT(I)  =  0.00039*YMV( I )  -14.5 

IF(HTT( I ).LT.0.0.  OR. DELT ( I ) . LT . 0 . 0 f  GO  TO  703 
IF(HTT(  n.GT.9.0.  OR.  DELT  ( I  )  .GT  .  9.0)  GO  TO  703 
CZ(I)  =  (  YMV(I)  -  YMV(I-l) 

HZd)  =  (  XMV(I)  -  XMV(I-1)3«*2 
NZ  =  NZ+1 

FA  =  2.*RA0IUS(I  )*SQRT(  OZ(I)  +  HZd)  ) 

SUM  =  SUM  +  FA 

IF(HTT( I ) .L7.2.125.AND.0ELT ( I ) . G7 . 2. 87 5 )  H7T( I )=2.125 
C  C0NVER7  70  THESIS  FORMAT  8.5X11. 

HTT(I)  =  hT7(I)*0.725 
OELT(  I  )  =  OELTd  )=f'0.725 

IF(  H7T  (I  )  .GE.5.6  .OR.  OELTd). GE.  6.5  )  GO  TO  4 
703  CONTINUE 

4  CALL  LINE{HTT,05LT,MZ,  1,  2  ) 

NM  =  NZ*KX 
DO  603  I=1,NM 

CCAA(I)  =  0.00039*CCA( I )  -  4. 

IF(  CCAA(I)  .LE.  0.0  )  CCAA(I)  =  0.0 
IF(  CCAA(I)  .GT.  7.3  )  CCAA(I)  =  7.5 
CCAA(I)  =  0.0C039*DDA(  I  )-  14.5 
IF(  CCAA(I)  .LE.  0.0  )  DOAA(I)  =  0.0 

I F(CCAA( I ) .LT.2. 12  5.AN0.0DAA( I ) .GT.2.875  )  CCAA( I  )  =  2.12 
IF(  CCAA(I)  .GE.  9.0  )  ODAA(I)  =  9.0 
C  CONVERT  TO  THESIS  FORMAT  8.5X11. 

CCAA(  I  )  =  CCAAd  )=«'0.725 
DCAA(  I)  =  DDAA( I ) *0.725 
603  CONTINUE 

SUM3  =  0.0001*SUM 
WRIT5( 6,304)  SUM3 
HRITE(6,303)  RADIUS(NZ),  NZ 
CALL  LINE(  XT,YT,KK,  1,  1  ) 

CALL  LIN£(CCAA,DOAA,NM,  1,-7  ) 

CALL  PL0T(0. 0,12. 0,-31 


666  00  10  1=1, KK 

XPd  )  =  .0306*XO(  I  ) 

YPd  )  =  .0306*Y0  (  I  ) 

10  CONTINUE 

DC  667  J=1,N 
AX(J)  =  AX(J)*0.78 
AY( J) =AY( J )*0.78 
667  CONTINUE 

00  668  K=1,M 
XR(K)  =  XR(K)*0.78 
YR(K)  =  YR(K)*0.78 
668  CONTINUE 

CALL  LINE(  XP,  YP,  KK  ,-1,1) 

CALL  SYMB0LC1.44, 10.0,.14, 

I'ANALYSIS  OF  OIL  SPILL  MOVEMENT ',  .0, 30  ) 

CALL  SYMeOL(1.75,  9. 7, .14, 'IN  SAN  FRANCISCO  BAY  AREA', 
2.0,25  1 

CALL  SYMB0L(0.77,  9. 4, .14, 

I'AS  PREDICTED  BY  ESTUARINE  NON  TIDAL  OR  I  FT '  ,0.0, 41  ) 

CALL  SYMBOL(5.8,  7 . 95 , 0 . 07 , ' SUI SUN  BAY', 0.0, 10) 

CALL  SYMBGL(0.20,  7.  74 , . 07 ,  ' 38N '  ,  0  .0 , 3 ) 

CALL  SYMBOL( 2.30, 7 .78, .07, ' SAN  OABLO  BAY', 0.0, 13) 

CALL  SYMBQL(5.06,  7. 54, .07,'*  BEN  IC I A '  , 0  .0 , 9 ) 

CALL  SYMB0L(3.70,5.60,.07,  '*  A L B ANY ' ,0 . 0 , 8  ) 

CALL  SYMBOL (4. 40, 4. 10  ,.14,  'ALAMEDA' ,0.0,7) 

CALL  SYMBOL(2.17,3.73,0.07,  'HUNTER  PC  I  NT  '  ,  0 . 0 , 1 2 ) 

CALL  SYMBOL{3.5,2.72,.07, 'SAN  FRANCISCO  BAY', 0.0,  17) 
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ooo 


CALL  SYMBOLO. 50, 1.75  ,.07, 'FOSTER  CITY  «',.0,13) 
CALL  SYMBOL (4,55 ,0.4,0.14, ' PALO  ALTO', 0.0, 9) 

CALL  L!N5(AX,AY,  N,  1,1) 

CALL  LINE(XR,YR, M,  1,1) 

CALL  PLOKO. 0,14. 0,-3) 

CALL  PLGTE 


STCP 
END 

//GO.SYSIN  DO  ’i' 


/ 
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A  P  P  E  N  D  I  X 


B 


COMPUTER  SIMULATION  PROGRAM  OF 
THE  SOCIAL  COST  OF  THE  LOSSES  IN  THE  RECREATION 
ACTIVITIES  AND  THE  COMMERCIAL  FISHING  ACTIVITIES 


//HAN1874  JOB  (2696,  1446, RL44)  ,  'iMUDJIARCJO  SMC  1874* 

//  EXEC  FORTCLGP, REGION. G0=180K 
//FGRT.SYSIN  DO  ’i' 

DIMENSION  XD(50),  YD(50),  P(50),  7X ( 10 ) , TY (  10 ) 
DIMENSION  XH(50>,  YH(50),  XC(50»,  YC ( 50  I ,  YCC(50) 
DIMENSION  XE(IO),  Y5(10),  PYY(IO),  V(10),  SC(IO) 

C  INITIAL  INPUT  DATA 

REA0(5,89)  ALFA, BEFA,PX,PY, AK, BK 
K  =  4 

C  READ  VOLUME  SPILLED  IN  LITERS. 

READ(5,83)  1V(J),  J=1,K) 

WRITE(6,85  ) 

WRITE(6,82)  (V(J),  J=1,K) 

C 

79  FORMAT!/////) 

80  FORMAT!/, 10X,'N0: (• ,13, •)• ,2X, 'SOCIAL  COST=',  F12.2, 
12X,'  VOL.  SPILLED  LTRS=',F12.1  ) 

81  FORMAT! / ,58X, • AH=' ,F9 .2, 3X, ' BH  =  '  ,  F9.2) 

32  FORMAT!  /,10X,5F12.1  ) 

83  FORMAT!  5F12,1  ) 

84  FORMAT!/, lOX, 'PYY!  ',  13,'  )  SO.  PRICE  LEVEL  OF  Y:',F8.2) 

85  FORMAT! /,4X, '  INPUT  DATA:  '  ) 

36  FORMAT! » 1' ,4X, 'COORD.  P~S  OF  PRGD.FNC,  MAX.  PROFIT', 

1'  LINE,  optimal  PROFIT  LINE  AND  WELFARE  FUNCTION'  ) 

87  FORMAT!//,  4X,'  EQUILIBRIUM  POINT:  '  ) 

88  FORMAT !5X, ' AL PHA  =  ' , F4. 2 , 3X , ' BET A= ' , F4. 2 , 3 X  ,  '  PX= '  , 
1F8.2,3X, 'PY=* ,F3.2,3X,' AK=' ,F9.2,3X, 'eK=' ,F9.2  ) 

89  FORMAT!  2F4. 2, 2F 8. 2, 2F9.2  ) 

90  FORMAT!/, 10X,F9.1,3X,F9.1,5X,F9*1,5X,F9.1,3X,F12.1  ) 

91  FORMAT!/, 5X, '$  PROFIT  OF  X  £  Y=  ',F12.2,3X, 

I'WELFARE  INITIAL  STATUS  :', F8. 1  ) 

92  FORMAT!  / , lOX , F 3 . 2 , 3 X , 4F9 . 1  ) 

93  FORMAT !/, lOX, 'X  UN  ITS  :  '  , F9. 1 , 3 X ,  '  Y  UNITS  :', F9. 1 , 3X , 
l'J=',I3  ) 

94  FORMAT!'!',’  COORDINATE  POINTS  OF  SUPPLY  AND  DEMAND', 
2'  CURVE'  ) 

N  =  25 

M  =  20 

CW  =  2.50 

XH!11  =  0.0 

ZETA  =  ALFA  +  BEFA 

DO  1004  J=1,K 

C  PRCDUCTION  FUNCTION  AS  A  FUNCTION  OF  VOLUME  SPILLED. 

AKl  =!V!  J)  /10.=<=*6  )  **2 

AH  =  AK*!0.9*£XP ! -AKl  )+  0.1) 

BKl  =  !5  .’i'V!  J  )/10  .*^6)=i'*2 

BH  =  BK*!0,75*EXP(-BK1)  +  0.25) 

WRITE!6,79) 

WRITE(6,88 )  ALFA, BEFA,  PX,PY,AK,BK 
WRITE(6,81)  AH,  BH 
XUP  =  PX*AH>i'=S‘2 
YUP  =  PY*BH=!‘*2 

QP  =  SQRT  !PX**2-AH-'!==i‘2  +  PY=<'=<'2’i=BH*=i'2) 

RL  =  XUP/QP 
RM  =  YUP/QP 
C  EQUILIBRIUM  POINT. 

XE(J)  =  SORT  (ALFA^i'AH^  +  Z/  ZETA  ) 

YE!J)  =  SQRT!BEFA*8H=^*2/  ZETA  ) 

C  MAXIMUM  PROFIT 

PBOFT=  PX=i'!XUP/OP )  +  PY*(YUP/gP) 
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C  WELFARE  STATUS  VALUE 

WWX  =  ALFA*PROFT/  (  PX=!=ZETA) 

WWY  =  BEFA*PROFT/  (  PY’i'ZETA) 

WS  =  CW^XE(  J  )*=!=ALFA  YE(J)**B£FA 
WRITE(6,87) 

WRITE(6,93)  RL,RM,  J 
WRITE(6,93)  WWX,WWY,  J 
WRITE(6,93)  XE(J),  YE(J),  J 
WRITE(6t91)  PPOFT,WS 

C  M3VING  OF  PRICE  TO  EQUILIBRIUM  STATE. 

6UL  =  BH**2  -  YE ( J 

IF(  BUL.LT.  0.0)  BUL=0.0 

PYY(  J)  =  YE(  J)*PX*AH/(  BH=i'SQRT(BUL  )  ) 

WRITE(6,84)  J,PYY(J) 

WRITc(6,86) 

DO  1000  1=1, N 

C  PPCDUCTICN  POSSIBILITIES  CURVE 

QRS  =  BH**2  -  (  BH^=^2/AH*’!'2  )*XH(  I  )=!'=«'2 
IF(  QRS  .LT.  0.0)  QRS  =  0.0 
YH4I)  =  SQRT(QRS) 

C  BUDGET  CCMSTRAINT  OR  PROFIT  FUNCTION 
YC(I)  =  (PROFT  -  PX*XH(i))/PY 
IF{  YC(I)  .LT.  0.0  )  YC(I)=0.0 

C  NEW  OPTIMAL  PROFIT  LINE 

PROFC  =  PX=!'XE  (  J  )  +  PYY(J)^YE(J) 

YCCd  )=(PROFC  -  PX=i'XH(  I  )  )/PYY(  J  ) 

IF(  YCCd  )  .LT. 0.0)  YCCd)=0.0 

C  WELFARE  CURVE 

Xhd  +  1)  =  XHd)  +  50. 

XOd  )  =  XHd  ) 

IF(  XD(I)  .EQ.  0.0  )  XD(I)  =  10.0 
YCd)  =  (WS/(CW^'XD(  I  l^^'^ALFA  )  )**<  1./3EFA  ) 

WRITE (6, 90  )  XH( I ) , YH( I ) ,YC ( I ) , YCC  (I ) ,YD(  I  ) 

1000  CONTINUE 
WRITE(6,94  ) 

C 

C  CONVERT  TO  THESIS  FORMAT 
DO  1001  I=1,N 

Xh(I)  =  0.00005*XHd  )  +  0.5 
IF(  XH( I )  .GT .6.5 )  XH(I)  =  6.5 
YFd)  =  0.00005=i‘YHd  )  +  0.5 
IF(  YH( I ) .GT.9.0)  YH(I)=9.0 
YCd)  =  0.0Q005=!'YC  (I  )  +  0.5 
IF(  YC(n. GT.9.0)  YC(I)=9.0 
XCd)  =  0.00005*XDd)  +  0.5 
IF(  X0(  I  )  .GT.6.5  )  XOd  )=6.5 
YDd  )  =  0.005=!'YDn  )  +  0.5 
IF(  YD( I ) .GT.9.0)  YD(I)=9.0 
YCCd)=  0.00005=i'YCCd  )+  0.5 
IF(  YCC  (  I  )  .GT.9.0  )  YCCd)  =  9.0 

1001  CONTINUE 
TXd)  =  0.5 
TY(1)  =  0.5 
TX(2)  =  6.5 
TY(2)  =  0.5 
TX(3)  =  6.5 
TY(3)  =  9.0 
TX(4)  =  0.5 
TY(4)  =  9.0 
TX(5)  =  0.5 
TY(5)  =  0.5 
CALL  PLOTS 

CALL  SYMBOLC  2 . 5 , 0. 0  ,  . 14, '  X  UNIT  PRODUCED  ',0.0,17  ) 
CALL  SYMBOL!  0 .0 , 2  .0  , . 14,  •  Y  UNIT  PRODUCED  ’,90.0,17  ) 
CALL  LINE!  TX,TY,5,1,1  ) 

CALL  SYMBOL!4.0,8.0,.07,  PROD.  PCSS  CURVE  ',0.,19) 
CALL  SYMB0L!4.0,7.3,.07,’-  MAX  PROFIT  LINE  ',0.0,18) 
CALL  SYMBOL! 4.0, 7.6, .07, WELFARE  CURVE  ',0.0,16) 

CALL  LINE!  XH,YH,  N,l,2  ) 

CALL  LINE!  XH,YC,  N,l,-7  ) 

CALL  LINE!  XH,YCC,N,1,6  ) 

CALL  LINE!  XD,YD,  N,l,5  ) 
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ooooo  o  oo  o  o  oo 


CALL  PLCT{  0.0, 12. 0,-3  ) 


P(l)  =  5. 

DC  1002  I=1,M 
SUPPLY  CURVE 

XH(I)  =  AH*=!=2*P(  I  ) /PROFC 
YHd)  =  6H*=i'2*P{  I  ) /PRCFC 
DEMAND  CURVE 

XCd)  =  ALFA«PR0FC/ZETA*d./P{  I  )  ) 

YCd)  =  B£FA=!'PR0FC/ZETA=i=(  l./P(  I  )  ) 

WRITE  (6 ,92)  P(  I  )  ,  XH(  I  )  ,YHd  )  ,XC  d  )  ,YCd  ) 

Pd  +  1)  =  Pd)  +  5. 

1002  CCNTINUE 

CONVERT  TO  THESIS  FORMAT 
DC  1003  I=1,M 
PCI )  =  0.025*Pd  )  +  0.5 
IF(  Pd  ).GT.9.0)  Pd  )=9.0 
Xhd)  =  0.00005*Xh(I  )  +  0.5 
IFCXHCI  )  .GT.6.5)  XHd  )  =  6.5 
XC(  I )  =  0.00005*XCd  )  +0.5 
IFCXCd  ). GT.6.5)  XCd  )  =  6.5 
YCd)  =  0.00005*YC  (I  )  +  0.5 
IFCYCd  ). GT.6.5)  YCd  )  =  6.5 
Yhd)  =  0.00C05*YHd)  +  0.5 
IF(YH( I ).GT.6.5>  YH(I )  =  6.5 

1003  CONTINUE 

CALL  SYMBOL!  2 .5  ,0.0  ,  . 14, *  X  UNIT  PRODUCED  ',0.0,17 
CALL  SYMBOL!  0 . 0 , 2 .0 , . 14, •  PRICE  PER  UNIT  PRODUCT  ' 
190.0,24  ) 

CALL  SYMBOL!4.0,3.2,.07, ' 

CALL  SYMB0L!4. 0,8.0,. 07, ' 

CALL  LINE!  XC , F ,  M,l,l  ) 

CALL  LINE!  XH,P,M,1,2  ) 

CALL  LINE!  TX,TY,5,1,1  ) 

CALL  PLOT!  0.0, 12. 0,-3  ) 


SUPPLY  CURVE  •  ,0.0, 16 ) 
DEMAND  CURVE  • ,0.0, 16 ) 


CALL  SYMBOL!  2 .5 ,0 .0, . 14, • 
CALL  SYMBOL!  0 .0 ,2 .0 , . 14, ' 
190.0,24  ) 

CALL  SYMBJL!4.0,8.2,.07, ' 
CALL  SYMB0L!4.0,8.0,.07, • 
CALL  LINE!  YC, P,  M,l, 1  ) 
CALL  LINE!  YH,P,  M,l,2  ) 
CALL  LINE!  TX,TY,5,1,1  ) 
CALL  PLOT!  0.0, 12. 0,-3  ) 
CALL  PLCTE 
1004  CCNTINUE 


Y  UNIT  PRODUCED  •  ,0.0, 17 
PRICE  PER  UNIT  PRODUCT  ' • 


+  SUPPLY  CURVE 
-  DEMAND  CURVE 


' ,0.0,16) 
•  ,0.0, 16  ) 


SOCIAL  COST  OF  RECREATION  ACTIVITIES  C  COMMERCIAL  FISHING 
ACTIVITIES  !AND  RELATED  INDUSTRIES) 

WRITE!6,79) 

DO  1005  J=1,K 

SC!J)  =  PX#!XE!1)  -  X£!J))  +  P YY ! 1 ) ^ ! YE ! 1 )  -  Y£!J)) 
WRIT£!6,80)  J,  SC!J),  V!J) 

1005  CCNTINUE 
STOP 
END 

//GO.SYSIN  DC  * 
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